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PREFACE 


The  DELFIC  (Defense  Land  Fallout  Interpretative  Code) 
Output  Processor  Module  has  been  updated  and  revised  to  process 
data  supplied  by  the  new  Diffusive  Transport  Module  (DASA  2669 
and  its  supplement) .  This  document  describes  the  revised  Output 
Processor  Module  code;  it  replaces  DASA-18G0-VI. 
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1.  INTRODUCTION 


This  document  is  intended  to  fulfill  two  needs:  (1)  to 

provide  information  to  the  person  who  is  interested  in  understanding 

the  Output  Processor  only  in  sufficient  detail  to  make  use  of  it, 

and  (2)  to  provide  a  detailed  explanation  of  the  Output  Processor  to 

the  researcher  or  programmer  who  would  make  modifications  or  additions. 

The  sections  entitled  "Program  Description"  and  "User  Information"  are 

intended  to  fulfill  the  first  need;  the  sections  "Program  Details"  and 

"FORTRAN  Listings,"  the  second  need. 

The  original  DELFIC  Outp>.  Processor  Module,  as  described 

in  DASA-1800-VI^ ,  is  designed  to  process  grounded  fallout  parcels 

(2) 

that  are  output  by  the  DELFIC  Transport  Module  .  These  parcels  are 
nuclear  cloud  subdivisions  in  the  form  of  square  wafers.  In  the  hori¬ 
zontal  plane,  they  possess  discrete  boundaries:  their  particulate 

content  is  uniformly  distributed  between  these  boundaries.  The  new 

(3) 

Diffusive  Transport  Module'  yields  radically  different  descriptions 
of  grounded  fallout.  Grounded  fallout  parcels  produced  by  the  Diffu¬ 
sive  Transport  Module  are  called  deposit  increments.  Each  deposit 
increment  is  distributed  in  the  ground  plane  via  a  bivariate  Gaussian 
function.  The  processing  requirements  for  the  deposit  increments  are 
sufficiently  different  from  those  for  the  square  wafers  that  most  of 
the  old  code  is  obsolete.  Therefore,  with  the  exception  of  a  few 
important  subroutines,  the  entire  code  has  been  rewritten. 

This  document  is  prepared  in  the  format  of  its  predecessor, 
and  many  parts  of  both  are  similar.  User  requirements  and  printed 
output  are  maintained,  where  possible,  in  their  original  fora.  The 
intended  applications  of  the  DELFIC  code,  which  are  to  provide  a 
numerical  research  tool  and  to  serve  as  a  fallout  prediction  standard, 
are  unchanged. 
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2.  PROGRAM  DESCRIPTION 

2.1  The  Purpose  and  Function  of  the  Output  Processor 

In  simplest  terms,  it  is  the  task  of  the  Output  Processor 
to  accept  descriptions  of  grounded  fallout,  process  the  deposit 
increment  data,  make  requests  for  particle  activities  or  mass  chain 

'  i.  \ 

concentral::.ons  from  the  Particle  Activity  Module1  when  required, 
accumulate  the  results  into  a  two-dimensional  memory  array  or  map 
image,  and  then  print  the  resulting  array  in  a  form  suitable  for 
viewing  as  a  map. 

The  code  provides  the  following  functional  capabilities 

1.  Great  flexibility  in  program  use  is  provided  in  terms  of  the 
variety  of  computations  that  are  available  (see  Table  4) . 

2.  The  Output  Processor  is  capable  of  handling  a  large,  essentially 
unlimited,  set  of  deposit  increment  data.  Since  this  set  can 
exceed  high  speed  access  memory  capacity,  an  open-ended  philoso¬ 
phy  is  adopted  for  its  treatment  by  use  of  peripheral  sforage 
capacity 

3.  The  area  coverage,  in  terms  of  location,  range,  and  scale  of  the 
map,  is  under  the  direct  control  of  the  researcher.  This  gives 
the  user  the  ability  to  produce  maps  for  superposition  on  other 
preexisting  maps:  and  it  enables  him  to  achieve  either  a  micro¬ 
scopic  or  a  macroscopic  view’  of  the  predicted  fallout  field. 

4.  The  Output  Proceosor  is  capable  of  handling  output  maps  that 
contain  a  larger  number  of  map  grid  points  than  can  be  stored 
in  the  computer  memory  at  one  time.  Thus,  the  code  has  open- 
ended  capability  with  regard  to  map  size. 

5.  In  computing  radiation  exposure  rates  at  arbitrarily  specified 

timjs,  it  is  deemed  of  great  importance  to  avoid  reliance  on  a 

-1.2 

single  time  decay  function  (such  as  t  ),  which  is  applicable 
only  to  a  mixture  of  unfractionated  fission  products  —  not  in 
general  to  isolated  samples  of  fallout  such  as  those  that  appear 
locally  in  fallout  fields.  Therefore,  the  Output  Processor 


►  vff  fffl  t  f 


Module  is  built  to  work  in  close  ’iaison  with  the  Particle  Acti- 

(4) 

vity  Module  so  that  activities  can  be  computed  directly  from 
the  primary  mass  chain  data  for  particles  at  the  particular 
time  or  times  specified  in  each  output  request.  Furthermore, 
and  as  a  consequence  of  this  approach,  the  user  can  request 
computation  and  display  of  concentrations  of  any  particular  mass 
chain. 

6.  With  regard  to  display  of  the  fallout  map  data  produced  by  the 
Output  Processor  we  are  faced  with  somewhat  conflicting  require¬ 
ments:  (1)  we  desire  a  numerical  display  of  the  data  rather 
than  some  sort  of  purely  pictorial  or  graphical  display  because 
of  the  intended  research  and  comparison  standard  applications  of 
the  system,  whereas  (2)  an  automated  pictorial  or  graphical  dis¬ 
play  relieves  the  user  of  the  txme  consuming  and  tedious  task 

of  coping  with  numerical  tabulations  and  hand  contouring.  The 
display  actually  provided  is  a  compromise.  A  numerical  display 
is  provided;  however,  it  is  in  a  format  that  allows  strips  of 
the  printed  computer  output  to  be  attached  side  by  side  so  that 
the  entire  fallout  prediction  area  is  included  on  the  assembled 
paper.  Thus,  the  output  tabulation  consists  of  the  requested 
output  data  printed  on  each  of  the  points  of  a  spatially  undis¬ 
torted  grid.  The  assembled  map  can  be  easily  contoured  directly 
on  the  printer  output  paper.  The  major  disadvantage  of  this 
type  of  display  is  that  sometimes  the  map  assemblies  are  quite 
large . 

7.  The  Output  Processor  is  simple  to  use  and  is  reasonably  foolproof 
ar.u  automatic  with  respect  to  its  internal  operations.  Since  the 
sizes  of  input  and  output  data  sets  can  vary  widely,  the  code 
contains  a  certain  amount  of  essentially  "dimension  free" 
programming. 
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2.2  Inputs  to  the  Output  Processor 

The  primary  input  to  the  Output  Processor  is  the  tape  of 
deposit  increment  descriptions  that  is  prepared  by  the  Diffusive 
Transport  Module  (see  Table  5) •  In  addition  to  the  deposit  incre¬ 
ment  descriptions,  this  tape  contains  Hollerith  identifiers  for  the 
preceding  DELFIC  module  runs,  and  a  collection  of  critical  data  such 
as  explosion  yield,  ground  zero  coordinates,  height  of  burst,  a  fall¬ 
out  particle  size  class  table,  etc.  The  data  set  for  each  deposit 
increment  consists  of  (see  Table  1) :  impact  coordinates  of  its 
center  of  mass,  its  impact  time,  its  particle  diameter  (each  deposit 
increment  is  composed  of  monodisperse  particles) ,  total  mass  of 
particles  in  it,  and  the  parameters  needed  by  the  bivariate  Gaussian 
function  to  distribute  it  in  the  ground  plane. 

In  addition  to  the  tape  input,  the  user  must  communicate 
to  the  program  via  card  input  his  wishes  regarding  types  of  output 
computations  and  map  specifications.  He  must  provide  run  identifi¬ 
cations.  And  he  must  supply  printer  characteristics  data  that  are 
necessary  for  production  of  undistorted  maps.  The  run  identifier  is 
an  arbitrary  72-character  Hollerith  statement  which  the  user  can  set 
to  identify  and  associate  outputs  and  inputs.  The  printer  charac¬ 
teristics  data  are  the  number  of  characters  per  inch  printed  by  the 
off-line  printer  in  the  cross-page  and  down-page  directions.  Map 
specifications  are  the  geographical  limits  of  the  map,  the  distances 
of  separations  between  map  points,  and  choice  of  format  for  printing 
individual  map  point  ordinate  values.  The  data  displayed  in  the  map 
are  for  one  of  of  the  options  listed  in  Table  4  and  discussed  below. 

2.3  Computation  and  Display  Options 

The  following  is  a  listing  and  brief  discussion  of  the 
major  options  for  computation  and  display.  An  exhaustive  list  of 
all  currently  available  options  is  provided  in  the  "User  Information" 
chapter.  (See  Table  4.) 


1.  Prirf^d  descriptions  of  impacted  particles 

Under  this  option  the  contents  of  the  deposit  increment  tape 
(IPOUT)  are  printed.  This  option  is  valuable  In  checking  the 
execution  of  experimental  transport  codes,  and  it  is  also  useful 
in  providing  a  hard  and  readable  copy  of  the  stored  results  of 
transport  production  runs. 

2.  Computation  options 

The  descriptions  below  apply  to  each  ordinate  value  of  a  map. 

a.  Count  of  contributing  deposit  increments.  This  can  be  of 
value  to  the  user  in  assessing  the  statistical  significance 
of  computed  quantities  at  all  points  cn  the  map. 

b.  Exposure  rate  ’'normalized"  to  time  H  +  1  hour*.  This  i-~ 
the  option  that  is  most  commonly  used  for  comparing  fallout 
patterns.  It  should  be  noted  that  differences  may  exist 
between  DELFIC  H  +  1  hour  normalizations  and  those  resulting 
directly,  or  indirectly,  from  backward  extrapolations  of 
field  data.  In  backward  extrapolations  a  single  decay 
constant  is  usually  used  through  the  map  area,  whereas 
DELFIC  provides  a  more  rigorous  modeling  of  radioactive 
decay . 

c.  Exposure  rate  at  time  II  +  Tl.  This  is  che  exposure  rate 
at  H  +  T1  taking  into  account  the  impact  times  of  all 
deposit  increments. 

d.  Exposure  accumulated  from  H  +  Tl  to  infinity.  This  is  the 
exposure  as  integrated  from  time  H  +  Tl  or  particle  impact 
time,  whichever  is  later. 


*  A  computation  of  radiation  exposure  or  other  quantity  that  is 

"normalized"  to  time  T  assumes  that  deposition  is  complete  through¬ 
out  the  map  area  at  time  T.  Thus  if  T  is  small,  the  normalized 
values  may  be  larger  than  actually  could  be  observed  at  that  time. 
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e.  Exposure  accumulated  from  time  H  +  Tl  to  time  H  +  T2.  This 
is  the  exposure  as  integrated  from  time  H  +  Tl  or  deposit 
increment  impact  time  (whichever  is  later)  to  time  H  +  T2. 

A  faster  alternative  treatment  of  accumulated  exposure  not 
accounting  for  deposit  increment  impact  time  is  also  provided. 

f.  Total  fallout  mass  per  unit  area  of  deposition  plane.  This 
is  the  mass  of  fallout,  both  radioactive  and  inert,  deposited 
on  the  map  grid  points  during  the  entire  fallout  period. 

g.  Fallout  mass  per  unit  area  deposited  between  times  Tl  and  T2. 
This  is  the  fallout  mass,  both  active  and  inert,  deposited 
during  the  specified  interval. 

2 

h.  Activity  produced  by  a  user  specified  mass  chain  (curies/m  ). 
Preparation  of  undistorted  maps 

The  Output  Processor  produces  a  numerical  presentation  of  fall¬ 
out  data  on  a  spatially  undistorted  grid.  The  user  must  supply 
map  grid  spacing  values  for  both  directions,  and  he  must  supply 
the  printer  characteristics  (characters/inch  both  cross-page 
and  down-page) .  The  program  automatically  adjusts  one  of  the 
grid  spacings  just  enough  to  accommodate  the  printer  character¬ 
istics  so  that  spatial  distortion  is  avoided.  A  map  produced  by 
the  Output  Processor  cc.~l.str  of  a  sequence  of  numbered  "strips" 
of  computer  printer  paper  which  can  be  assembled  side-by-side 
into  a  single  map  of  the  overall  area  covered.  When  so  assembled 
the  data  point  with  minimum  x  and  minimum  y  coordinates  will  be 
found  in  the  lower  left-hand  corner  of  the  map  (i.e.  the  lower 
left-hand  corner  of  strip  number  one) .  The  coordinates  of  this 
point  will  be  (XMIN  +■  DGX,  YMIN  +  DGY) .  This  point  need  not  be 
either  the  origin  of  coordinates  or  ground  zero. 

Numerical  display  formats 

Two  options  exist  at  this  time  for  printing  ordinate  values  at 
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the  map  grid  points.  These  options,  which  we  designate 
as  the  two-line  E  format  and  the  two-line  F  11.3  format,  are 
explained  and  illustrated  as  follows  for  a  single  map  ordinate: 


NNNNNN 

+  V.VW, 

which  is  to  be  interpreted  as 
±  NNNNNNV.VW. 


2.4  Deposit  Increment  Processing 

The  Output  Processor  prepares  all  maps  with  their  grid 
points  aligned  in  the  west-east  and  south-north  directions.  The  map 
x  coordinate  direction  is  positive  toward  the  east,  and  the  y  direc¬ 
tion  is  positive  toward  the  north.  The  z  coordinate  direction  is 
positive  upward.  Each  deposit  increment  is  defined  by  the  data 
listed  in  Table  1. 

The  standard  deviations,  a,  ,  c,,  and  angle,  a,  are  calc»  - 

lated  by  the  Diffusive  Transport  Module  they  are  unique  for 

(3)  2 

each  fallout  parcel  .  c„  is  the  variance  of  the  Gaussian  distri¬ 
buted  deposit  increment  in  the  average  downwind  direction.  It  is  the 
sum  of  the  initial  value  input  from  the  Cloud  Rise-Transport  Interface 

Module ^  and  the  downwind  component  of  the  turbulent  dispersion 

2 

variance  that  is  computed  for  the  parcel  trajectory.  oJl  is  the 
corresponding  crosswind  variance.  The  angle  a  is  the  angle  between 
the  positive  x  axis  and  the  average  downwind  direction  axis.  The 
averaging  is  a  space-weighted  averaging  computed  along  the  trajectory. 
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TABLE  1 

DEPOSIT  INCREMENT  DESCRIPTION  PARAMETERS 


Mathematical 

symbols 

FORTRAN 

Mnemonics 

Parameter  Definition 

x  ,  y  ,  z 

F  7P  P 

X(I),  Y(I), 
ZOUT(I) 

space  coordinates  of  the  center  of 
mass  (meters)* 

t 

P 

T(I) 

time  of  deposit  (seconds) 

0M  ,  °1 

SXOT(I) , 
SYOT(I) 

Gaussian  distribution  standard  devia¬ 
tions  in  the  (horizontal)  downwind 
and  crosswind  directions  (meters)* 

a 

ROUT(I) 

angle  between  the  downwind  direction 
axis  and  the  positive  x  axis 
(radians)* 

D 

PS(I) 

fallout  particle  diameter  (micrometers) 

M 

FMAS (I) 

mass  of  fallout  (kilograms; 

See  text  for  a  -sore  complete  definition. 


The  vertical  coordinate  of  a  deposit  increment,  z  , 

P 

usually  is  the  same  as  that  of  the  deposition  plane.  However,  for 
fallout  parcels  that  are  advected  through  one  of  the  vertical  wind- 
field  boundaries,  or  that  are  not  impacted  when  the  transport  time 
boundary  is  reached,  the  OTM  records  the  z^  at  the  level  of  boundary 
penetration.  The  output  processor  code  rejects  any  deposit  incre¬ 
ment  whose  Zp  is  ten  meters  or  greater  above  the  deposition  plane. 

Consider  a  deposit  increment  with  total  mess  or  activity 
content  Q.  Then  at  a  point  x,y,  the  areal  density  of  mass  or  acti¬ 
vity  q(x,y) ,  is 


where 


q(x,y) 


_  Q- 

2*o, ,  o± 


exp 


(1) 


X  =  x  cosa  +  y  sina  (2) 

Y  =  y  cosa  -  x  sina  (3) 


and  Xp  and  are  defined  similarly.  The  %.y,  and  X,Y  .oordinate 

axes  are  related  as  shown  in  Figure  1. 

An  input  datum  to  the  code  is  a  parameter  QCUT  =  qm--n 

that  represents  a  threshold  value  for  a;  I  deposit  increments.  At 

any  map  point,  a  contribution  from  any  deposit  increment  that  is  less 

than  q_ .  is  ignored.  To  provide  efficient  processing  of  deposit 
Tu.n 

increments,  we  need  a  simple  and  fast  method  for  determining  the 
boundary  that  encloses  q(x,y)  >  q^  for  individual  deposit  incre¬ 
ments.  The  method  used  is  discussed  next. 

Let  q(x,y)  equal  q^  and  take  logarithms  of  both  sides 
of  Eq.  (1).  Then  we  get 


For  y  >  0,  Eq.  (4)  is  an  equation  of  an  ellipse.  This  ellipse  is 
drawn  in  Figure  1.  We  wish  to  determine  its  tangent  lines  parallel 
to  the  y  and  x  axes,  which  are  labeled  x^.x^,  and  y^yj,  in  the  figure. 
From  the  general  properties  of  an  ellipse,  it  can  be  shown  that 
these  lines  are 


V*T 
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2  2  2 
cos  a  +  ov  sin"  a). 


and 


yT.y^. 
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2  2  2  2 
2y  (o„  sin  a  +  o,  cos  a) 


(6) 


(7; 


In  che  code,  computation  of  x^,x^  and  y^y.^  *s  use<*  to  establish 
whether  or  not  a  deposit  increment  contributes  to  a  particular  map 
or  map  section. 

For  each  deposit  increment,  map  points  are  considered 
row-by-row.  fne  bounding  rows  aie  determined  from  the  values. 

For  a  particular  row,  the  bounding  x  coordinates,  xc,x^  are  given  by 


.  (8) 


3E 


2.5  Processing  for  Impact  Time  and  Particle  Size 

Maps  can  be  prepared  for  time  of  fallout  onset,  time  of 
cessation,  smallest  particle  deposited  and  largest  particle  deposited. 
For  these  options,  the  deposit  increment  processing  methods  described 
in  the  preceding  section  are  applied  as  follows.  At  a  specified  map 
point,  a  particular  deposit  increment  is  considered  or  bypassed 
depending  on  whether  or  not  the  map  point  lies  within  the  contribu¬ 
tion  ellipse  defined  by  Eq.  (4).  In  computing  y  (Eq.  (5)),  for 

mass  per  unit  area  is  used.  When  the  map  point  falls  within  the 
contribution  ellipse,  the  impact  time  or  particle  size  of  the  deposit 
increment  may  be  rejected  or  it  may  replace  the  value  already  stored 
for  the  map  point,  depending  on  the  outcome  of  a  straightforward 
logical  test. 

2.6  Sequences  of  Processing  Requests 

The  Output  Processor  accepts  in  a  single  input  a  sequence 
requests  for  processing.  The  user  can  obtain  any  number  of  maps 
any  descriptions  in  a  single  run  provided  that  the  same  deposit 
increment  input  tape  (IPOUT)  is  used  for  them  all.  The  code  is 
completely  open-ended  in  this  respect. 

The  utility  of  this  feature  is  illustrated  by  the  follow¬ 
ing  example.  Suppose  the  user  desires  various  maps  to  be  prepared 
for  each  of  two  different  sets  of  map  specifications.  For  example, 
the  user  may  desire  large-scale  maps  of  essentially  the  entire  local 
fallout  field  for  (1)  exposure  rate  normalized  to  H  +  1  hour, 

(2)  total  accumulated  exposure,  and  (3)  activity  from  mass  chain  95. 

Es  may  also  desire  these  options,  plus  some  o^srs,  for  a  h-’gb  reso¬ 
lution  map  that  covers  a  geographically  smaller  area  close-in  to 
ground  zero.  Tc  accomplish  this  he  can  specify  the  map  limits  and 
grid  intervals  for  the  large-scale  map  and  follow  it  by  the  needed 
computation  option  request  cards.  These  data  would  be  followed  in 
turn  by  the  other  nap  specifications  and  another  series  of  computation 
option  request  cards. 
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2.7  Output  Processing  Independent  of  Other  DELFIC  Modules 

In  its  primary  role  the  Output  Processor  ects  as  the 
terminal  module  of  the  DELFIC  system.  Nevertheless,  it  can  also 
operate  independently  of  the  other  programs  of  the  DELFIC  system 
except  for  the  Particle  Activity  Module.  This  feature  can  be  used 
to  advantage  if  the  user  saves  the  magnetic  tape  results  of  the  trans 
port  program's  execution.  Thus,  the  user  need  not  specify  all 
desired  output  at  the  time  of  the  transport  execution  but  can  make 
subsequent  runs  of  the  Output  Processor  as  specific  questions  arise 
during  the  course  of  his  research.  The  tape  and  card  :  iputs  to  the 
Output  Processor  are  the  same,  regardless  of  which  way  the  program 
is  used. 

2.8  General  Logic  of  the  Output  Processor 

In  this  section  we  present  a  cursory  description  of  the 
operations  of  the  Output  Processor,  including  organizational  flow 
charts.  More  thorough  descriptions,  which  include  detailed  discus¬ 
sions  of  the  more  involved  subroutines,  are  given  in  the  Program 
Details  chapter. 

The  Output  Processor  subroutines  are  listed  with  brief 
functional  descriptions  in  Table  2.  In  addition  to  these  programs, 
a  contro.l  program  and  the  utility  subroutine  ERROR  are  required.  A 
control  program,  OPP,  which  was  used  for  independent  operation  of 
the  Output  Processor  Module  on  the  UNIVAC  1108  computer,  is  included 
in  die  FORTRAN  listings.  The  FORTRAN  listings  of  OPP  and  ERROR  are 
self  explanatory. 

Output  Processor  operations  are  separated  into  two  main 
parts;  these  are  controlled  by  subroutines  LINK8  and  LINK9.  LINK8 
(Figure  2)  is  used  for  run  initialization.  It  also  can  be  used 
solely  to  print  the  contents  of  a  DTM  binary  output  tape,  IPOUT.  If 
maps  are  to  be  created,  LINK8  calls  the  initialization  portion  of  the 
Particle  Activity  Module  (PAM)  code,  PAM1.  PAMl  prints  out  the  PAM 


TABLE  2 


OUTPUT  PROCESSOR  PROGRAM  SYNOPSIS 


Program 

Name 


Purpose 


LINK8 


LINK9 


PC HECK 


Initializes  and  writes  printout  headings.  Prints  con¬ 
tents  of  tape  IPOUT  if  requested.  Calls  first  part  of 
Particle  Activity  Module  (PAM1)  to  perform  request- 
invarient  part  of  activity  calculations. 

Controls  request-dependent  portion  of  the  Output 
Processor  computations.  Calls  the  second  part  of 
the  Particle  Activity  Module  (PAM2). 

Accumulates  contributions  from  individual  deposit 
increments  into  the  map  point  ordinates. 

Controls  flow  of  deposit  increment  description  data 
blocks  to  and  from  tape. 

Prints  the  fallout  maps. 

Initializes  for  a  map  calculation.  Cbmputes  deposit 
increment  contribution  boundaries  in  the  map. 


PDMP 


Sorts  out  deposit  increments  tha.  will  contribute  to 
subsequent  map  core  loads  (if  any)  and  dumps  them  onto 
tape  for  temporary  storage. 


*  **tr*&* «  **s*  „  ^ 


START 


MAKE  TAPE  ASSIGNMENTS 
AND  SET  ARRAY  COUNTERS 


READ  HEADER  DATA  PROM  THE 
BINARY  INPUT  TAPE,  IPOUT 


READ  RUN  IDENTIFICATION  \ 
AND  CONTROL  ARRAY  1C  ] 


PASS  TAPE 
IPOUT 


/P  rintX 
/deposit  \ 

/INCREMENT^ 

.DESCRIPTIONS 


PASS  TAPE  IPOUT  1 
AND  PRINT  THE  DEPOSIT 
INCREMENT  DESCRIPTIONS . 


'  READ  PRINTER  N 
SPECIFICATIONS  AMO 
PRINT  OUTPUT 
,  HEADER  j 


/ areX 

/  THE  \ 

/  DEPOSIT  \ 
/INCREMENTS  \ 
JO  BE  PROCESSED 
\  INTO  MAPS  / 


STOP 


CALL  PAM1 


RETURN 


RETURN  TO  THE 
DELFir  MONITOR 
FOR  A  CALL  TO 
LINK9 


FIGURE  2.  ORGANIZATIONAL  FLOW  CHART  CF  SiBROUTINE  LINKS 


'-  i&* ‘jf^355*^£^ 


t 


i 


« 


htiader,  and  performs  the  shot-specific  but  map-independent  portion 
of  the  PAM  computations. 

LINK9  (Figure  3)  accepts  map  specifications  and  map 
requests;  it  controls  the  processing  of  the  deposit  increment  data 
into  maps  and  the  printing  of  the  maps.  Deposit  increment  data  are 
read  from  the  DTM  binary  output  tape,  IPOUT.  A  complete  pass  of 
tape  IPOUT  is  made  for  each  map  that  is  created.  The  map-specific 
portion  of  the  Particle  Activity  Module,  PAM2,  is  called  for  maps 
that  require  activity  calculations. 

A  map  specification  defines  map  boundaries  and  grid 
intervals.  An  unlimited  number  of  map  specifications  can  be  accom¬ 
modated.  For  each  map  specification  an  unlimited  number  of  map 
requests  can  be  accommodated.  A  map  request  selects  one  of  the  six¬ 
teen  computation  options  that  are  available  (see  Table  4) ,  and 
provides  data  that  are  specific  for  that  request.  With  reference  to 
Table  3  ,  map  specifications  are  input  via  cards  4,  5,  and  6,  and 
map  requests  via  cards  7. 

Core  storage  of  map  ordinate  data  is  carried  in  the 
singly  dimensioned  array  OMAP  (see  card  124  in  the  LINK8  listing) . 
The  dimension  of  OMAP  must  correspond  to  the  value  assigned  to  the 
variable  NMAP  (see  card  182  in  the  LINK8  listing) .  When  the  OMAP 
array  is  not  large  enough  to  accommodate  an  entire  map,  the  program 
will  still  function,  provided  that  two  scratch  tapes,  JPOUT  and 
KPOUT  (see  cards  178  and  179  of  the  LINK8  listing)  are  provided. 

By  use  of  these  scratch  tapes,  maps  with  essentially  unlimited 
numbers  of  points  can  be  prepared. 


- 


FIGURE  3: 


ORGANIZATIONAL  FLOW  CHART  OF  SUBROUTINE  LINK9 
AND  ASSOCIATED  PROGRAMS 
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3.  PROGRAM  DETAILS 

In  this  chapter  we  discuss  each  of  the  subroutines  listed 
in  Table  2  in  some  detail.  These  discussions,  in  conjunction  with 
the  organizational  flow  charts  (Figures  2  and  3) ,  are  intended  to 
provide  the  background  needed  to  easily  grasp  the  complete  logical 
and  computational  content  of  the  FORTRAN  code. 

3.1  Subroutine  LINK8 

This  subroutine  initializes  for  an  OPM  run.  One  pass 
through  LINK8  is  made  for  each  run.  LINR8  can  be  used  to  initialize 
for  map  preparation,  or  it  can  be  used  simply  to  print  the  contents 
of  the  (binary) deposit  increment  tape,  IPOUT,  that  has  been  prepared 
by  the  Diffusive  Transport  Module  (DTM) .  If  maps  are  to  be  prepared, 
LINK8  calls  subroutine  PAM1  of  the  Particle  Activity  Module  (PAM)  v 
PAMl  performs  the  portions  of  the  particle  activity  calculations  that 
are  not  specific  to  individual  map  requests. 

The  operations  of  LINK8  are  particularly  simple.  Figure  2 
presents  an  essentially  complete  outline  of  them. 

3 . 2  Subroutine  LINK9 

Subroutine  LIKK9  initializes  for  individual  map  specifi¬ 
cations  and  map  requests.  It  begins  by  reading  in  a  map  specification 
data  card.  These  data  are  the  map  boundary  coordinates,  the  ^p  grid 
intervals,  and  a  combined  ground  roughness  and  instrument  response 
factor  for  gamma  radiation  exposure  rate.  If  the  sum  of  absolute 
values  of  the  grid  intervals  is  found  to  be  zero,  the  run  is  terminated. 
If  not,  the  altitude  relative  to  mean  sea  level  of  the  fallout  deposi¬ 
tion  plane  is  read  in.  Then  an  integer  control  array,  JC,  is  read  in. 
This  is  used  to  specify  the  map  ordinate  numerical  display  format. 
Finally,  these  input  data  are  printed. 

A  map  request  card  is  read  and  the  OMAP  array  is  initialized. 
If  the  computation  option  code  parameter  NRF.Q  (see  Table  4)  is  zero. 
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the  end  of  map  requests  for  this  map  specification  is  signaled;  con¬ 
trol  is  transferred  to  read-in  of  the  next  map  specification  card. 

If  NREQ  is  greater  than  zero,  the  map  request  data  are  printed.  If 
an  activity  map  is  requested,  some  additional  initialization  tasks 
are  performed  and  PAM2  is  called.  PAM2  computes  the  particle  acti¬ 
vity  array  FP,  which  contains  the  total  activity,  in  suitable  units 
according  to  the  request,  associated  with  each  particle  size  class. 

Next,  the  map  grid  intervals  are  adjusted,  if  necessary, 
to  provide  an  undistorted  printed  map.  The  printer  characterization 
parameters  IV  and  IH  are  used  to  make  this  adjustment.  Then,  the 
number  of  map  points  in  the  x  and  y  directions  for  the  complete  map 
are  computed.  If  the  whole  map  cannot  be  accommodated  in  core 
storage,  the  number  of  core  loads  beyond  the  first,  NZ,  is  computed. 

If  NZ  =  0,  subroutine  GOGO  is  called  to  begin  processing 
the  deposit  increment  data  into  the  map.  When  all  of  the  deposit 
increment  data  have  been  processed,  control  is  returned  to  LINK9  and 
subroutine  MAP  is  called  to  print  the  map. 

If  NZ  >  0,  the  map  must  be  prepared  in  NZ  +  i  parts*. 

Two  scratch  tapes,  KTAPE  and  LTAPE,  are  used  for  temporary  storage  of 
deposit  increment  description  data  that  contribute  to  succeeding  map 
parts*.  These  tapes  are  treated  like  the  DTM  output  tape,  IPOUT, 
when  they  are  used  as  input  for  successive  map  part  computations. 

For  NZ  >  0,  the  code  begins  by  initializing  the  KTAPE  and 
LTAPE  assignments.  Then  it  calls  GOGO  and  MAP  to  prepare  and  print 
the  first  part  of  the  map.  Next,  it  enters  a  loop  indexed  from  1  to 
NZ  in  which  the  remaining  map  parts  are  prepared  and  printed. 


These  map  parts  should  not  be  confused  with  the  printed  map  strips 
that,  when  assembled,  constitute  a  complete  map.  What  we  have 
called  a  map  part  constitutes  a  portion  of  a  map  that  can  be  con¬ 
tained  in  the  computer's  rapid  access  memory.  In  general,  each 
such  map  part  will  yield  more  than  one  map  strip. 
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Finally,  the  sun  of  all  map  ordinates,  FSUM,  is  printed, 
and  control  is  passed  to  read-in  of  the  next  map  request  card. 

3.3  S  ub  r ou  t ine  GOGO 

Subroutine  GOGO  reads  into  core  storage  a  block  of  deposit 
increment  description  data  from  binary  input  tape  KTAPE.  Tape  KTAPE 
is  either  the  DTM  output  tape,  IPOUT,  or  one  of  the  two  scratch  tapes 
used  for  temporary  storage  when  a  complete  map  cannot  be  contained  in 
core. 

Each  data  block  is  preceded  on  the  input  tape  by  an  integer 
block  count  NIJ.  When  NIJ  =  0,  this  signals  that  the  end  of  the  input 
tape  has  been  reached.  For  NIJ  >  0,  the  block  of  deposit  increment 
data  is  read  into  core.  Subroutine  PCHECK  is  then  called  by  GOGO  to 
process  the  block  of  deposit  increment  data  into  the  map. 

On  return  of  control  to  GOGO,  the  values  of  parameters  NZ 
and  ICTR  are  compared.  NZ  is  the  number  of  map  core  loais  (map  parts*) 
beyond  the  first  required  to  prepare  the  map.  Both  NZ  and  ICTR  are 
set  by  subroutine  LINK9.  If  NZ  -  ICTR,  no  additions'  map  core  loads 
are  required,  and  another  block  of  data,  preceded  by  its  block  count, 
is  read  in  from  tape  KTAPE.  If  NZ  4  ICTR,  a  succeeding  map  core  load 
is  signaled.  In  this  case,  subroutine  PDMP  is  called.  PDMP  writes 
on  to  temporary  storage  tape  those  deposit  increment  descriptions 
that  will  contribute  to  subsequent  map  parts.  Then,  the  next  block 
of  data  is  read  from  tape  KTAPE. 

3. A  Subroutine  PDMP 

When  an  entire  map  cannot  be  contained  in  core  storage  in 
the  OMAP  array,  the  map  must  be  prepared  in  two  or  more  parts  via 
construction  of  two  or  more  map  core  loads*.  In  this  case. 


*  See  footnote,  page  20. 
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after  each  block  of  deposit  increment  descriptions  is  processed  into 
the  current  in-core  map  part,  subroutine  PDMP  is  called  by  subroutine 
GOGG  to  write  on  to  temporary  storage  tape  the  data  for  those  deposit 
increments  that  will  contribute  to  subsequent  map  parts. 

During  the  processing  of  the  deposit  increments  into  the 
in-core  map  part,  subroutine  PCHECK  labels  each  deposit  increment  to 
indicate  whether  or  not  it  will  contribute  to  subsequent  map  parts. 
This  labeling  is  done  in  array  RTR  (see  the  PCHECK  glossary  in  the 
FORTRAN  listings).  PCHECK  also  tallies  the  number,  NE,  or  deposit 
increments  currently  stored  in  core  that  do  not  contribute  to  subse¬ 
quent  map  parts. 

The  first  operation  in  PDMP  is  to  compute  the  number,  KP, 
of  deposit  increment  descriptions  that  must  be  saved.  Next,  the 
storage  block  of  depos.'t  increment  descriptions  is  rearranged  so  that 
all  of  the  data  to  be  sived  are  stored  in  a  continuous  block  in  the 
low-core  end  of  the  arrays.  Finally  the  block  count  followed  by  the 
block  of  deposit  increment  descriptions  are  copied  out  onto  tape 
LTAPE. 

3.5  Subroutine  PCHECK 

Subroutine  PCHECK  is  called  by  subroutine  GOGO  to  initiate 
the  processing  of  a  core-stored  block  of  deposit  increment  data  into 
the  core-stored  map  or  map  part. 

The  subroutine  operations  are  wholly  enclosed  in  a  DO 
loop  that  passes  the  complete  block  of  stored  deposit  increment 
descriptions.  The  discussion  that  follows  applies  to  each  deposit 
increment  in  the  block. 

First  the  altitude  of  the  deposit  increment  is  compared 
with  that  of  the  fallout  deposition  plane.  If  the  deposit  increment 
is  ten  meters  or  more  above  the  deposition  plane,  it  is  rejected  for 
further  processing.  If  not,  processing  continues. 

On  the  basis  of  the  value  of  NREQ,  the  computation  option 
code  (see  Table  4 ) ,  control  is  transferred  to  an  appropriate  portion 
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of  code  to  initialize  for  the  deposit  increment  processing.  This 
initialization  established  the  value  of  F,  which  is  equivalent  to  Q 
in  Eqs.  (1)  and  (5). 

Below  statement  number  100,  the  deposit  increment  contri¬ 
bution  ellipse  boundaries  (see  Figure  1),  XPRMU,  XPRML  and  YPRMU, 

YPRML,  are  computed.  These  are  equivalent  to  the  x^,  x^,  and  y^,  y^, 
values  given  by  Eqs.  (6)  and  (7).  The  ellipse  boundaries  are  tested 
against  the  map  boundaries  to  determine  if  the  deposit  increment 
contributes  to  the  map.  The  ellipse  boundaries  also  are  checked  to 
determine  if  the  deposit  increment  will  contribute  to  subsequent  map 
parts,  if  any.  The  deposit  increment  is  labeled  accordingly  via 
array  KTR.  If  the  deposit  increment  contributes  to  the  currently 
stored  map  part,  subroutine  CALC  is  called  for  further  processing. 

If  not,  further  processing  is  bypassed. 

3.6  Subroutine  CALC 

Subroutine  CALC  is  called  by  subroutine  PCHECK  to  compute 
contributions  from  an  individual  deposit  increment  to  the  map  ordinates 
and  add  or  enter  them  into  the  OMAP  array. 

The  first  operation  is  the  computation  of  those  factors 
and  terms  in  Eq.  (8)  that  are  independent  of  individual  map  point 
coordinates.  Next,  the  bounding  map  row  indices  NOB  and  NOT,  are 
computed.  These  are  established  by  YPRMU  and  YPRKL  values  (y^  and 
y p  that  have  been  computed  by  subroutine  PCHECK.  The  remainder  of 
the  operations  are  contained  in  a  loop  that  is  indexed  between  NOB 
and  NOT. 

On  each  map  row,  as  defined  by  its  row  index,  all  points 
have  a  common  y  coordinate.  Therefore,  in  each  pass  through  the  row 
index  loop,  the  limiting  x  coordinates,  x^,  x^,  in  the  row  are  deter¬ 
mined  by  application  of  Eq.  (8).  On  the  basis  of  the  xc,  x^,  the 
limiting  map  colmun  indices,  NOL  and  NOR,  in  the  row  are  computed. 

Then,  the  OMAP  array  index  extremes  for  points  in  the  row,  K  and  L, 
are  computed,  and  an  inner  loop  indexed  between  K  and  L  is  entered. 
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Within  this  inner  loop,  contributions  to  the  OMAP  array  elements  are 
computed  and  added  to  the  elements  or  replace  them  depending  on  the 
requirements  of  the  computation  option. 

The  OMAP  array  is  singly  dimensioned.  The  map  points  are 
represented  in  the  array  in  the  following  order.  The  array  element 
0MAP(1)  represents  the  lower  left-hand  corner  point  in  the  map.  The 
array  is  then  filled  by  the  successive  points  in  the  lowermost  row. 
Following  che  rightmost  point  in  the  lowest  rov?  in  the  leftmost  point 
in  the  next  to  lowest  row,  and  so  on. 

3.7  Subroutine  MAP 

This  subroutine  writes  map  print  images  on  the  operating 
system  output  tape,  ISOUT.  It  writes  a  map  title,  a  description  of 
the  quantity  that  the  map  portrays,  and  an  indication  of  the  ordinate 
format  used.  It  divides  rhe  output  map  into  printer  strips  on  the 
basis  of  the  parameter  INC,  which  is  the  number  of  map  ordinate 
columns  that  can  be  accommodated  by  the  printer  paper.  It  prints  a 
strip  count  (MAPRUN)  at  the  top  of  each  strip  for  identification 
purposes.  A  separate  call  of  MAP  is  necessary  for  each  map  part  or 
core  load*. 

Following  the  FORTRAN  statement  listing  of  subroutine  MAP, 
we  see  at  its  beginning  a  transfer  to  a  first-pass  portion  of  code 
if  MAPRUN  equals  zero.  In  this  first-pass  portion  of  code,  parameter 
initializations  are  performed,  a  map  title  is  written,  the  display 
option  control  parameter,  JC(1),  is  checked  for  an  acceptable  value, 
and  then  a  brr.ich  transfer  is  made  to  a  code  that  writes  the  ordinate 
format  identification  and  makes  control  transfer  assignments  for  use 
within  the  map  writing  loops. 

Between  the  statement  numbers  102  and  170  a  two-part 
title  is  written  that  describes  the  quantity  presented  in  the  map. 
Between  statement  numbers  170  and  2023  initializations  are  made  for 
the  three  nested  map  writing  loops.  When  2023  is  first  reached,  M 

*  See  footnote  page  20. 
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contains  the  number  of  printer  strips  that  are  to  be  produced,  and 
LEFT  has  the  number  of  columns  that  should  appear  on  the  last  printer 
strip. 

At  2023,  which  is  the  return  point  for  the  outer  map  writ¬ 
ing  loop  (printer  strip  loop),  MAPRUK,  the  counter  of  printer  strips, 
is  incremented  and  the  strip  title  is  written.  Also,  KL,  the  lower 
index  for  retrieval  from  the  one-dimensional  map  array  OMAP,  is  set 
at  its  initial  value.  Note  that  in  the  iteration  KL  progresses  from 
its  largest  value  to  its  smallest  value  to  invert  the  map  which  is 
stored  numerically  inverted  in  the  map  array. 

At  card  number  229,  the  return  point  for  the  middle  map 
writing  loop  (printer  line  loop),  KH,  the  upper  index  for  retrieval 
from  the  map  array,  is  set  and  KDC,  an  index  for  the  printer  line 
integer  array  JMAP,  is  initialized. 

At  card  number  237,  the  return  point  for  the  last  map 
writing  loop  (data  point  loop) ,  KDC  Is  incremented  and  a  transfer  is 
made  to  the  desired  presentation  format  code  on  the  basis  of  previous 
assignment.  The  two  printer  format  codes  take  their  inputs  from  the 
map  array  and  place  their  results  back  into  the  map  array  and  into 
the  integer  printer  line  array  JMAP.  All  map-producing  codes  return 
to  statement  number  300. 

Below  300  the  printer  lines  are  written  onto  the  output 
tape,  certain  indexing  operations  are  performed,  and  return  is  made 
to  deal  with  either  the  next  line  in  the  current  strip  or  the  first 
line  (and  title)  on  the  next  strip,  or  a  final  return  is  made  to  the 
calling  program.  Note  that  if  entrance  is  made  to  MAP  with  MAP RUN 
set  positive  as  a  consequence  of  a  previous  entrance,  the  overall 
titles  will  not  be  printed  again  and  strip  counting  will  be  resumed 
where  it  had  been  left  off. 
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4.t  USER  INFORMATION 

This  chapter  is  intended  to  be  useful  as  a  user’s  manual 
fjr  the  Output  Processor  Module  (0PM).  However,  a  reading  of  this 
chapter  alone  is  not  sufficient  preparation  for  use  of  the  0PM;  the 
user  also  should  at  least  read  sections  2.1  through  2.3. 

Inputs  to  the  0PM  are  of  two  kinds:  (1)  a  card  input 
that  contains  the  user's  specifications  for  individual  fallout  maps, 
and  (2)  a  binary  tape  input  that  contains  fallout  deposit  increment 
descriptions  prepared  by  the  Diffusive  Transport  Module  (DTM)  as  well 
as  other  information  passed  on  by  the  preceding  DELFIC  modules. 

4.1  Card  Input 

Card  inputs  are  listed  in  loading  sequence  in  Table  3. 

Notice  that  the  Particle  Activity  Module  (PAM)  card  Inputs  are 

(4) 

imbedded  in  this  deck.  The  user  is  referred  to  DASA-1800-V  for 
descriptions  of  these  cards. 

As  indicated  in  Table  3  ,  the  map  requests  are  segre¬ 
gated  into  sets.  Each  map  request  set  is  introduced  by  a  map 
specification  card  (card  4)  that  defines  the  map  boundaries  and  grid 
intervals.  Each  map  request  card  (card  7)  that  follows,  specifies 
a  map  option  (Table  4 )  for  a  particular  map  that  is  to  be  computed 
and  printed.  Each  request  set  is  terminated  by  a  blank  card.  The 
run  is  terminated  by  a  blank  map  specification  card.  Thus,  the  last 
two  cards  in  the  deck  must  be  blank. 

CARD  1  -  Run  Identification 

A  description  of  the  0PM  run  is  input  via  this  card. 

CARD  2  -  Run  Control  Variables 

This  card  allows  up  to  eighteen  control  variables  to 
be  input.  Currently,  only  two  of  these  are  used: 
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TABLE  3 


CARD  INPUTS  FOR  IDENTIFICATION 


AND 

CONTROL  OF  THE  OUTPUT  PROCESSOR 

Card 

FORTRAN  Mnemonic 

Data  Set 

No. 

Content 

and  format 

1 

Run  identification 

OPID(J) ,  J=l,12 
(12A6) 

Initialization 

2 

Run  control  variables 

IC(J),  J=l,18 
(1814) 

3 

Printer  characteristics: 

IH,IV  (214) 

Number  of  characters  per  inch 
in  the  cross-page  and  down- 
page  directions 


Particle  Acti¬ 
vity  Module 
Card  Deck 


4 

Map  specification  data: 

XMAX,  XMIN, 

maximum  and  minimum  x  coor- 

YMAX,  YMIN, 

dinates,  maximum  and  minimum 

DGX,  DGY, 

y  coordinates,  grid  intervals 
in  the  x  and  y  directions 

GRUFF  (7F10.3) 

First  Set 

(all  in  meters) ,  a  combined 

of  Map 

ground  roughness  and  radia- 

Specification 

Cards 

tion  meter  response  factor 

5 

Deposition  plane  altitude 
(meters  relative  to  mean 
sea  level) 

ZDEP  (F10.3) 

6 

Map  control  variablas 

JC(J),  J=l,18 
(1814) 

(continued  on  next  page) 
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Table  3  (continued) 


Card  FORTRAN  Mnemonic 

Data  No.  Content  and  format 


7a 

Map  request:  computation  option 

NREQ,  Tl,  T2, 

code,  NREQ  (see  Table  A),times 

MASCHN,  QCUT, 

of  onset  and  cessation  of  the 

CUTMAP  (14, 

computation  (hr)  or  particle 

2F10.3, 

diameter  range  limits  (micro- 

14,  2F10.3) 

meters) ,  mass  chain  number , 

First  Set 

deposit  increment  contribution 

of  Map 

threshold,  map  ordinate  threshold. 

Request 

Cards 

7b 

Map  request 

7c 

Map  request 

• 

• 

7n 

• 

Map  request 

8 

Request  termination  blank  card 

Second  Set 

4' 

Map  specification  data 

of  Map 

Specification 

5’ 

Deposition  plane  altitude 

Cards 

6’ 

Map  control  variables 

7a' 

Map  request 

Second  Set 

7b’ 

Map  request 

of  Map 

• 

Request 

• 

Cards 

• 

8' 

Request  termination  blank  card 

Additional  Sets 

• 

of  Map  Specifi- 

• 

cations  and 

• 

requests 

• 

9 


Run  termination  blank  card 


Computation 
Code  NREQ 

0 

1 

2 

3 

4 


6 

7 


10 

11 

12 

13 

14 

15 

16 


TABLE  4 

HAP  COMPUTATION  OPTIONS 

_ Computation  Option  Description _ 

Termination  of  map  request  set 

Count  of  deposit  increments  contributing  to  each  map 
ordinate 

Exposure  rate  normalized*  to  time  H  +  1  hour 

Exposure  rate  at  time  H  +  T1  hours 

Integrated  exposure,  H  +  T1  to  im.j...  ..y  accounting 
for  time  of  arrival 

Integrated  exposure,  H  +  T1  to  H  +  T2  accounting  for 
time  of  arrival 

Total  mass  per  unit  area 

Total  mass  per  unit  area  deposited  from  time  H  +  T1 
to  H  +  T2 

Integrated  exposure,  H  +  T1  to  H  +  T2  assuming  all 
particles  have  arrived  by  H  +  T1  hours 

Same  as  8  integrated  to  infinity 

Activity  per  unit  area  from  an  individual  mass 
chain  (curie s/m^) 

Time  of  onset  of  fallout 

Time  of  cessation  of  fallout 

Smallest  particle  size  deposited 

Largest  particle  size  deposited 

Mass  per  unit  area  from  particles  in  size  range  T1 
to  T2 

H  +  1  hour  normalized*  exposure  rate  resulting  from 
particles  in  size  range  T1  to  T2  microns 


*  In  a  calculation  normalized  to  time  H  +  T,  it  is  assumed  that  ell 
fallout  is  grounded  at  time  H  +  T,  regardless  of  whether  this 
actually  is  the  case. 
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IC(17)  Controls  the  processing  of  deposit  increments. 
IC(17)  >  0  causes  the  program  to  stop  without 
entering  PAM1  or  LINKS).  This  setting  is  used 
if  only  a  printing  of  the  deposit  increment  tape 
(IPOUT)  is  desired,  or  if  the  user  simply  wants 
to  see  the  run  identifiers  for  the  preceding 
DELFIC  modules.  IC(17)  =  0  causes  a  normal 
entrance  to  the  main  body  of  the  Out-put  Processor 
regardless  of  whether  the  deposit  increment  tape 
has  been  printed. 

IC(18)  Controls  the  option  to  print  the  complete  con¬ 
tents  of  the  deposit  increment  tape,  IPOUT. 

IC(18)  >  0  causes  the  deposit  increment  tape  to 
be  printed.  IC(18)  -  0  bypasses  the  printing 
of  the  deposit  increment  tape.  For  either  zero 
or  positive  value  of  IC(18),  preceding  DELFIC 
module  run  identifiers  and  other  vital  run 
statistics  are  printed. 

CARD  3  -  Printer  Characteristics 

To  prepare  and  print  spatially  undistorted  maps,  the  Output 
Processor  needs  constants  which  describe  the  character 
spacing  of  the  off-line  printer  to  be  used.  These  con¬ 
stants  IH  and  IV  give  respectively  the  cross-page  and 
down-page  character  spacings  of  the  printer  in  characters 
per  inch.  If  IH  and  IV  are  found  to  be  zero,  the  program 
assigns  the  standard  values  of  10  and  6  to  them. 

CARD  4  -  Map  Specification  Data 

Maps  must  be  completely  specified  by  the  user.  He 
must  specify  limiting  coordinates  and  grid  intervals 
(grid  point  spacing).  All  maps  are  rectangular  in  shape 
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and  north-south,  east-west  in  orientation,  with  north 
always  at  the  top.  The  variables  XMAX  and  XMIN  indi¬ 
cate  respectively  the  maximum  and  minimum  values  of 
the  east-west  coordinates  of  the  map.  The  positive  x 
direction  points  eastward,  which  is  cross-page  to  the 
right  on  the  map.  YMAX  and  YMIN  similarly  indicate 
maximum  and  minimum  values  of  the  north-south  map 
coordinates.  The  positive  y  axis  points  northward, 
which  is  up-page  on  the  map. 


The  variables  DGX  and  DGY  indicate  the  map  grid-point 
separations  in  the  east-west  and  north-south  directions, 
respectively.  It  should  be  noted  that  on  the  printed 
map  the  actual  physical  spacing  of  the  data  points  is 
determined  in  part  by  the  printer's  character  and  line 
spacings.  Thus,  if  necessary,  the  code  uses  the  printer 
description  parameters,  IH  and  IV,  to  adjust  DGX  or  DGY 
so  that  a  truly  undistorted  map  is  produced.  In  per¬ 
forming  this  adjustment  the  program  uses  either  DGX  or 
DGY  as  the  scale  factor  basis,  depending  upon  which  of 
these  two  parameters  will  yield  the  largest  undistorted 
map  (smallest  scale  factor). 


Gamma  ray  exposure  or  exposure  rate  is  computed  for  a 
detector  placed  three  feet  above  an  unbounded  plane 
source  of  fallout.  The  computed  values  are  corrected 
for  ground  roughness  absorption  and  radiation  meter 
response  via  multiplication  by  the  factor  GRUFF.  GRUFF 
is  the  product  of  the  ground  roughness  attenuation 
factor  with  the  radiation  meter  response  factor.  A 
value  of  0.5  is  satisfactory  for  most  work. 


»iJrs. 
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CARD  5  -  Deposition  Plane  Altitude  | 

The  value  of  ZDEP  should  be  the  same  as  ZMIN  used  in  | 

the  Diffusive  Transport  Module  calculations.  j 

CARD  6  -  Map  Control  V  -.riables 

Currently  only  one  of  these  variables  is  used.  It  is 
used  to  specify  the  printing  format  of  the  map  ordinate 
values  (see  pp.  7  and  8). 

JC(1)  =  1  results  in  printing  of  the  map  ordinates  with 
a  two-line  E  format,  which  has  the  power  of  ten  printed 

on  one  line  and  the  associated  multiplier  printed  i 

immediately  below.  ! 

JC(1)  =  2  results  in  the  printing  with  a  two-line  F11.3 
format,  which  has  the  six  highest  order  characters 
printed  on  the  first  line  and  the  five  lowest  order 
characters  on  the  second  line. 

CARD  7  -  Map  Request 

The  computation  option  codes  (NREQ  values)  are  given 
in  Table  4  .  Except  for  options  15  and  16,  T1  and 
T2  represent  time  limits  (hours)  for  the  calcula¬ 
tions.  T1  is  the  earlier  time.  The  T1  field  may  be 

left  blank  for  NREQ  =1,  2,  6,  and  10  -  14.  The  T2  , 

field  may  be  left  blank  for  NREQ  =1-4,  6,  9,  and  ! 

10  -  i».  | 

| 

For  options  I?  and  16,  T1  and  T2  (T1  <  T2)  represent  i 

the  particle  diameter  extremes  (micrometers)  for  a  i 

range  of  particle  size  classes.  j 

j 

If  option  10  is  selected,  a  value  for  MASCHN  must  be  , 

specified.  It  is  the  atomic  number  of  the  radioactive  j 

mass  chain  for  which  output  is  to  be  displayed.  ] 


\ 
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The  quantity  OfUT  corresponds  to  q^  in  Eq.  (5).  It 
represents  the  threshold  value  of  the  area  density  of  the 
quantity  to  be  displayed  below  which  contributions  from 
individual  deposit  increments  are  to  be  neglected. 

The  quantity  CUTMAP  is  a  map  ordinate  threshold.  After 
all  contributions  have  been  accumulated  at  each  map 
point,  a  pass  is  made  through  the  map  ordinate  array 
(OMAP)  and  any  ordinate  with  a  value  less  than  CUTMAP 
is  set  to  zero. 

Values  for  QCUT  and  CUTMAP  that  have  been  found  adequate 
are: 


NREQ 

qcuT 

CUTMAP 

-4 

-2 

2 

10  4 

10 

-6 

-3 

6 

10 

10 

4 . 2  Tape  Input 

The  Output  Processor  requires  a  binary  tape  input  (unit 
IPOUT)  t..at  is  prepared  by  tne  Diffusive  Transport  Module.  This 
tape  contains  critical  run  data  and  run  identifiers  for  each  of  the 
preceding  DELFTO  modules.  It  also  contains  deposit  increment  des¬ 
criptions,  which  are  the  major  o.'tput  of  tha  Diffusive  Transport 
Module.  The  contents  of  tape  IFOUT  are  described  in  Table  5. 

4.2  Output 

An  example  of  the  OPK  output  is  given  in  the  -'’Sample  Printout" 
chapter.  The  Particle  Activity  Module  output  would  appear  between 
pp.  73  and  74  below. 

Note  on  p.  7 A  he  two  columns  of  numbers.  These  are  y  axis 
coordinates  that  are  printed  on  tha  same  scale  as  the  map.  They  can 
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be  cut  from  the  page  and  attached  along  the  sides  of  the  map  to 
specify  its  y  axis  coordinate  values. 

Units  of  quantities  displayed  in  the  maps  are:  ■ 
exposure  -  roentgens 
exposure  rate  -  roentgens  per  hour 
mass  per  unit  area  -  kilograms  per  square  meter 
time  -  seconds 

particle  diameters  -  micrometers 

activity  per  unit  area  -  curies  per  square  meter 


TABLE  5 


BINARY  TAPE  INPUT  TO  THE  OUTPUT  PROCESSOR 


Record 

No. 


Content 


Varinble  Names 


1  Tape  identification  word,  IPOUT 

2  Fission  yield  (KT),mass  of  the  cloud  soil 
burden  (kg),  soil  solidification  tempera¬ 
ture  (°K) ,  time  at  which  the  cloud  reached 
the  soil  solidification  temperature  (sec), 
geometric  standard  deviation  of  the  log¬ 
normal  particle  diameter  volume-frequency 
distribution,  total  yield  (KT) ,  altitude 
of  burst  above  mtUm),  x  coordinate  (E-W) 
of  GZ(m),  y  coordinate  (N-S)  of  GZ(m), 
detonation  time  (sec),  spare  data  word, 
fallout  particle  density  (kg/nr*),  the 
horizontal  cloud  subdivision  parameter 
IRAD  (see  Reference  5,  p.  56  ff.),  maxi¬ 
mum  cloud  radius  (m) ,  altitude  of  ground 
zero  above  msl(m). 

3  Cloud  Rise-Transport  Interface  Module 
run  identification 

4  Cloud  Rise  Module  run  identification 

5  Initial  Conditions  Module  run  identification 

6  Diffusive  Transport  Module  run  identification 

7  Number  of  particle  size  classes 

8  Particle  size  class  tables:  central 
particle  diameter  (pm),  volume  (mass) 
fraction,  particle  diameter  at  the  upper 
boundary  of  the  size  class  (pm) 

9  Number  of  (altitude)  entries  in  the 
atmosphere  description  tables 

10  Atmosphere  tables:  Altitude  relative  to 
msl(ra),  viscosity  (kg/m-sec) ,  density 
(kg /nr) 


JPOTJ 

FW,  SSAM,  SLDTMP, 
TMSD,  SIGMA,  TW, 
H3URST,  XGZ,  YGZ, 
TGZ,  BZ,  ROPART, 
IRAD,  RADMAX, 
ZBRSTZ 


FSEIU(I),  1=1, 12 

CRID(I) ,  1=1,12 
DETID(I),  1=1,12 
WID(I),  1=1,12 
ITaB 

PSIZE(I),  FMASS(I) 
PACT  ( I)  ,  1=1,  ITAB 

NAT  (=256) 


ALT I I) ,  AIEMP(I), 
RHO(I) ,  1=1,  NAT 


Table  5  (continued) 

Record 

Wo. 


11 

12 


13 

14 

«► 

15 


Content 


Deposit  increment  description  block  count 

Block  of  deposit  increment  descriptions. 
For  each  deposit  increment:  x,  y,  z,  and 
t  coordinates  (m  and  sec),  horizontal 
downwind  and  crosswind  dispersion  standard 
deviations  (m),  average  ’ind  heading  from 
due  east  (radians),  siz<-  class  central 
diameter  (pm),  mass  of  fallout  (kg) 

Block  count 

Block  of  deposit  increment  descriptions 


Zero  block  count 


Variable  Names 


NIJ 

X(I),  Y(I) ,  ZOUT(I), 
T (I) ,  SXOT(I) . 
SYOT(I) ,  ROUT(I), 

PS (I),  FMAS(I), 

1=1,  NIJ 


NIJ=0 
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5.  FORTRAN  STATEMENT  LISTINGS 


Except  for  the  control  program,  OPP,  which  is  given 
first,  the  subroutines  are  arranged  in  alphabetical  order  accord¬ 
ing  to  their  names.  A  listing  of  the  utility  program  ERROR  is 
given  in  DASA-2669(3) . 
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OPR 


OUTPUT  PROCESSOR  MAIN 
COMMON  /SET1/  nf411f„.. 

1CAY  »OETID(i2 )  ,OIAM«201I 

2FMASSC200I »IOISTR  ,IFXEC 
1NOSTR  , PS  f  2  0  0 )  »sr 

4T<iP2  ,T2M  ,US0IL 

57SCL  t NHOOO  ,ZV<200> 

OTMENSION  NUMT API 15  > 


,QMEAN 
,IRISE 
,SSAM 
,VPR 


, VX1200T 


,ONS 
,ISIN 
»T  <E 

t* 

,VY<200> 


,EXPO 
,ISOUT 
,TMPi 
, HEIGHT 


ISIN=S 
ISOUT =6 
TRISE=10 
DOl  I=1»15 
NU^TAP (I) =0 
NUMT  AP 1 2) =13 
NUMTAP<3)=1A 
NUMTAp(t,|  =  Q 
NUMTAP(9)=10 
NUMTAP(li)=ll 
NUMT  AP (12)  =  12 
CALL  LIN<3(NUMTAP» 

LIN*9 
EXIT 


CALL 

CALL 

STOP 

END 


OPP 
0  PP 
,OPP 
,OPP 
,OPP 
,OPP 
OPP 
OPP 
0°p 
OPP 
OOP 
OnP 
0°P 
OPP 
OPP 
OPP 
OPP 
OPP 
OPP 
OPP 
OPP 
OPP 
OOP 
OPP 


1 

2 

■* 


4 

5 

6 

7 

8 
9 

10 

11 

1? 

1? 

14 

19 

IE 

17 

18 

19 

20 
21 
22 

23 

24 
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SUBROUTINE  calc 

CALC 

1 

c 

CALC 

? 

c 

MA»A  VERSION  DFSICNFC  TO  OPERATE  WITH  THE  DTM 

CALC 

3 

c 

H.G.NORMEN'T  JUNE  25,19X1 

calc 

4 

c 

calc 

S 

r 

thts  suproutinf  computes  hap  contributions  FOR  INOIVIDUAL 

CALC 

6 

c 

PflRTIC 

LFS 

calc 

7 

c 

CALC 

e 

r 

CALC 

Q 

c** 

**««V**** 

**calc 

10 

c 

calp 

n 

r 

NOB 

SMALLEST  FOSSTOLE  y  INOFX  OF  a  CONTRIBUTION  FROM  A 

CALC 

l? 

r 

“APTIOLE 

Ffl  Lr 

13 

c 

NOl 

SMALLEcT  F0SSI9LE  X  INOFX  rp  A  CONTRIBUTION  F“OM  A 

CALC 

14 

c 

RAPTICIE 

calc 

is 

c 

N0° 

LARGEST  POSSIBLE  X  INOFX  0“  A  CONTRIBUTION  FROM  A 

calc 

i  ; 

r 

“article 

CALC 

17 

r. 

NOT 

LART.rci  PCSS  I0Lr  Y  INOFX  OF  A  CONTRIBUTION  FROM  A 

calc 

18 

r 

“ARTICLE 

calc 

10 

r. 

Y“EL 

y  roORflHATE  OF  THE  MAP  POINT  ROW  CUPRFNTLY  BE1“S 

calc 

20 

o 

C0NSTpFRrr  “ELATTVF  to  thf  parttcle  Y  COOPOTNS t  ■■ 

cal'' 

21 

n 

X“FL 

X  COORflNATE  OF  The  MAP  POINT  CUP»ENTLY  BEING 

calf 

Z2 

c 

C OnS T OFRf r  RELATIVE  TO  THF  PARTTCLE  X  COOPOINATF 

calc 

23 

c 

XL 

LEFT  a OUNO°Y  X  COOROINATE  OF  THE  PARTICLE 

CALC 

24 

c 

CONTRTFUTICN  ELLIPSE  IN  THE  YREL  MAR  “OW 

calc 

?c 

c 

X° 

RIGHT OCUNCRY  x  COORDINATE  OF  THE  “ARTICLE 

calc 

?f 

c 

C  ON T“  T FC T  ICN  ELlIRStr  jN  THE  YREL  HA °  ROW 

calc 

?7 

n 

NWX 

NUmofR  CF  MAR  PMTNTS  SOflNNc0  BY  A  PARTICLE 

calf 

?? 

r 

rONCFNTRATTON  ELLIPSE  IN  A  RCW 

CALC 

f  - 

c 

v  a  “  x  2 

2.o*oalsstan  OTSTbn.  vapiance  along  a  axis 

TALC 

3_r- 

c 

VARY? 

R.O’GALSSIAN  oistbn.  VARIANCE  along  b  axis 

CALC 

c 

C 

•1AGNITLOE  (I.E.  INTEGRATED  VALUE*  OF  a  PARTICLE 

CALC 

32 

c 

property  TO  “E  niSTRIBUTEO  CN  THE  MAP 

calc 

3x 

c 

calc 

?  4 

c 

AL^o  S 

Er  L I N< 3  GLOSSARY  ANO  PCHEC<  GLOSSARv 

CALC 

7C 

c 

calc 

36 

c  * 

**  calc 

37 

r 

calc 

3  P 

COMMON 

/Sc  T 1/ 

calc 

?B 

icay 

,  CrTTD(12l  .CIAM1201*  ,0MFAN  ,DNS  .EXPO 

.calc 

4  0 

2DITTO(?00)  ,niSTR  .TFXEC  .IRISE  ,ISIN  .ISOUT 

.calc 

41 

TNn^tP 

,TTO<200  »  ,sr  . S  S AM  .THE  ,TMP1 

,caic 

42 

e»TMP? 

»T?M  ,0  » VP“  ,W  .HBURST 

.calf 

•♦3 

SSCLOHP 

.  KHQOO  , ZV (200  *  , V  X (2 0  0 *  .VY (200* 

calc 

44 

COMMON 

/PARRAT/ 

calc. 

4C 

IX (500) 

,Y(500)  . 7CUT (500 *  ,SXOT(S00l  ,SYOT(500)  ,ROUT(500» 

.calc 

46 

2PS(500>  , FmAS (SOI)  ,<TO(SnO)  ,  c  ,  GAMA  ,  BSO 

.CALF 

47 

3&SO 

>  SIN A  . COSA  .WFMftS  (POO) .YBRWU  » YPRML 

.CALC 

43 

4T {SOOT 

CALC 

49 

COMMON 

7RUN0AT7 

CALF 

sc 

INI  J 

,NE  . NPEO  . N  Z  , IC  TP  .NXMAP 

.calc 

51 

2T1 

,T?  .MAPRUN  ,T-Z  . IP  , JC  (13) 

.calc 

52 

3IC(18> 

.NYMAP  .NTASK  ,N0“0  ,XGZ  ,YG7 

CALC 

53 

Common 

/maptat/ 

CALC 

54 

lOHAPil 

5000)  .OCUT  .CUTMAP  ,OGX  .OGY  .DELTAx 

.CALC 

55 

2  XMA  y 

,XMIN  .  YMAX  .  YMIN  ,FSUM  .R'JFSAH 

.CALC 

56 

3X1 

» X2  .M6TAPE  .ZnEP 

calf 

57 

<-» 

calc 
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c 
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c 

c 


c 
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c 

c 

c 
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c 

c 

c 


c 

c 
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100 

110 

120 

130 

140 


150 

160 


oara  prpgrn/6hcslc  / 

INITIALIZE  FO»  THIS  PARTICLE 

VAPX2=  ASQ/GAMA 
V ARY2=  PSO/GAMA 

A  =  SINA*COSA*(1.0/VAPY2-  1.  0 / VARX 21* 2. 0 
P  =  4. Q/VAPX2/VAPY2 

C  =  (CCSA**2/VAPX2  ♦  SINA*»2/VAPY2»*2. 0 
0  =  2.0*GAMA*C 

0=F/SX0T{IP>/SY0Tt IF) '6. 28 31 8531 

COMPUTE  SMALLEST  Y  INDEX  OF  A  CONTRIBUTION 

NOP  =  (YPRML  -  YMIN1/PGY 
NOt,  =  NOB+-l 

IFCNOP.LT. 11  N09=l 
IF  CNOB.LE. NYMaP)  GO  TO  120 
IRROR=-110 
GO  TO  400 

COMPUTE  LARGEST  Y  INDEX  OF  A  CONTRIBUTION 
NOT  =  {YPRMU  -  YMINJ/flGY 

ifcnot.gt.nymapi  not=nymap 
IFCNOT.GT.O  )  GO  TO  140 
IPROR=-130 
GO  TO  400 

ENTER  THE  MAP  ROW  LCCP 
00  350  J=NOB, NOT 

COMPUTE  THE  LIMITING  X  COORDINATES  OF  THE  PARTICLE  CONTRIBUTION 
ELLIPSE  IN  THIS  ROW 

YREL  =  J 

YRFL  =  YMIN  «•  DGY*YREL  -  YtIP» 

RADIC  =  -9*YREL**2*C 
IFCRA0IC.GE.0.01  GO  TO  160 
IRPOR=150 
RADIC=0. 0 

CALL  ERROR CPROGRH, IRRORfISOUTl 
RAOIC=SORT  (RADIO 
XL  =  XCIPJ  «•  (YREL* A-  RAOTC) /C 
XR  =  XL  ♦  2.0*RADIC7C 

COMPUTE  SMALLEST  X  INOEX  OF  A  CONTRIBUTION 

NOL  =  (XL-X11/DGX 
N0L=N0L*1 

IF(N0L.LT.1>  NOL=l 
IF(NOL.GT.NXMAP)  GO  TO  350 

COMPUTE  LARGEST  X  INOEX  OF  A  CONTRIBUTION 
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CALC 
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66 

calc 
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68 

calc 

69 

CALC 

7  0 

CALC 

71 

CALC 

72 

CALC 

73 

CALC 

74 

CALC 

75 

CALC 

76 

calc 

77 

CALC 

7? 

CALr 

79 

calc 

90 

calc 

81 

CALC 

82 

CALC 

8? 

CALC 

84 

calc 

95 

CALC 

36 

CALC 

87 

CALC 

88 

CALC 

89 

CALC 

90 

CALC 

91 

CALC 

92 

CALC 

93 

CALC 

94 

CALC 

95 

CALC 

96 

CALC 

97 

calc 

98 

CALC 

99 

CALC 

100 

CALC 

101 

CALC 

102 

CALC 

103 

CALC 

104 

CALC 

105 

CALC 

106 

CALC 

107 

CALC 

108 

CALC 

109 

calc 

110 

CALC 

111 

CALC 

112 

CALC 

113 

CALC 

114 

CALC 

115 

CALC 

116 

42 


non  non 


160  NOR  =  (XR-XU/OGX 

IF(NOR.GT.NXMAP)  N0R=NXMAP 
IF(NOR.LT.l)  GO  TO  350 
200  NWX  =  NOR  -  NOL  f  1 

IF(NWX,GT.0»  GO  TO  220 
IF(NWX.fQ.O)  GO  TO  250 
210  IRROR=-210 
GO  TO  <*00 

COMPUTE  OMAP(M)  ARRAY  INOFX  EXTREMES  FO®  MAP  POINTS  IN  THIS  ROW 

220  MC°MT  =(J-11*NXMAP 
<  =  NOL  4-  MCRMT 

L  =  XfNWX-1 

AO JUST  OR  A OO  CONiR  I8UTI0NS  TO  THE  MAP  POINTS 

GO  TO  (224, 22 4,221,221, 222, 222)  ,NORC 

221  OMA  =T ( IP) 

GO  TO  22V 

222  OMA -PS (IP) 

224  OO  300  M=K , L 

GO  TO  (225,245,230, 240, 230*240) ,NORD 

225  0MAP(H)=0MAP(M)4-1.  0 
GO  TO  300 

230  OMAPCM1  =  AMIN1 (OMA  ,CMAP(M) ) 

GO  TO  300 

240  OMflP(M)  =  AMAX1  ( OM  A  ,Gi'i AP(M)  I 
GO  TO  300 

245  XREL=I1  -MCRMT 

XREL  =  XI  4-  OGX*  XFEL  -  X(IP) 

OMA  =  Q*EXP(  -  (XREl*COSA  4-  YREL*S INAl **2/VARX2  -  (YREL»COSA 
1  -  XPEL*SINA|**2/VARY21 
250  OMAP ( M 1  =  OMAP(M)  4-  OMA 
300  CONTINUE 
350  CONTINUE 
RETURN 

400  CALL  ERRORtPROGRM, IFROR,ISOUT1 
END 
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CALC  127 
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CALC  133 
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CALC  136 
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CALC  142 
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SUBROUTINE  GOGO 

THIS  SUBROUTINE,  WNICH  IS  CALLED  BY  LINK9,  CONTROLS  REAO-IN  OF 
PARTICLE  DATA  ANO  ROUTS  ITS  PROCESSING  ANO  ITS  LOADING  ON  TO 
TEMPORARY  STORAGE  TAPE. 


H.G.NORMENT  JUNE  28,1971 


GOGO  8 


GLOSSARY 


GOGO 

10 

ICTR  A  CONTROL  PARAHETEP  -  WHEN  ICTR.NE.NZ  , 

ANOTHER 

GOGO 

11 

MAP  CORE  LCAO  IS  SIGNALED  TO 

FOLLOW 

GOGO 

12 

NIJ  A  BLOCK  CCUNT  OF  DATA  STORED 

ON  TAPF  AND/OR  IN  CORE 

GOGO 

13 

N7  NUMBER  OF  MAP  CORE 

LOAOS  REQUIRED  BEYOND 

THE  FIRST 

GOGO 

14 

GOGO 

15 

ALSO  SEE  LIN<9  GLOSSARY 

GOGO 

16 

GOGO 

17 

**GOGO 

18 

GOGO 

19 

COMMON  /SFT 1/ 

GOGO 

20 

1CAY  »DETID(12)  , DIAMT201) 

, DMEAN 

,ONS 

»EXPO 

,  GOGO 

21 

2BITIO(20Q) ,IOISTR  ,IEXEC 

, IRISE 

,  IS  IN 

, ISOUT 

,  GOGO 

22 

3N0STP  , TIOI20D)  , SD 

,  SSAM 

,TME 

,TMP1 

,  GOGO 

23 

4TMP2  , T2M  ,U 

»  VPR 

,H 

, HBURST 

,  GOGO 

24 

5SCLDHB  , NHODO  ,ZV(200) 

,VXt*200) 

»  VY (200) 

GOGO 

25 

common  /parcat/ 

GOGO 

26 

1XC500I  ,Y(500I  , ZOUT  (5  0  0) 

,  SXOT  (5001 

,SYOTI50C) 

,ROUTt500l 

,  GOGO 

27 

2PST500)  ,FNAS(500>  ,KTR<500) 

♦  F 

,  GAMA 

,  BSQ 

,  GOGO 

28 

3ASQ  , SINA  , COSA 

,HFMAS (200) ,YPRHU 

,  Y°RML 

,  GOGO 

29 

4T(500> 

GOGO 

30 

COMMON  /RUNOAT/ 

GOGO 

31 

INI J  , NE  , NPEQ 

»NZ 

,  ICTR 

,NXHAP 

,  GOGO 

32 

2T1  ,T2  , MAPRUN 

,  TGZ 

»IP 

, JC (18 1 

,  GOGO 

33 

3IC(18)  ,NYMAP  , NT ASK 

,NORD 

,XGZ 

,YGZ 

GOGO 

3^ 

COMMON  /CONOAT/ 

GOGO 

35 

1IP0UT  , JPOUT  , KPOUT 

,KTAPE 

,LT  APE 

,MARRAY 

,  GOGO 

36 

2NMAP  , MXREO  , IH 

,IV 

GOGO 

37 

GOGO 

39 

**GOGO 

39 

GOGO 

40 

OATA  PR0GRH/6HG0G0  / 

GOGO 

41 

IEXEC  =  1 

GOGO 

42 

REAP  A  DATA  BLOCK  CCUNT 

GOGO 

43 

GOGO 

44 

100  RE AO  CKTAPEI NIJ 

GOGO 

45 

NE  =  0 

GOGO 

46 

GOCO 

47 

ARE  HE  FINISHED  PROCESSING  THE 

DATA- 

GOGO 

49 

GOGO 

49 

IF(NIJ.EO.O)  GO  TO  400 

GOGO 

50 

IF (NIJ.LE. M ARRAY)  GC  TO  200 

GOGO 

51 

150  IRROR=-150 

GOGO 

52 

160  CALL  ERROR (PROGRM, IRROR,ISOUT) 

GOGO 

53 

GOGO 

54 

REAO  A  BLOCK  OF  PARTICLE  BATA 

GOGO 

55 

GOGO 

56 

200  REAO(KTAPE) (Y( II , Y ( I) ,7CUT II) , T Cl) ,SXOT(I) 

, SYOTC I) ,ROUT (I) »PSfII 

,  GOGO 

57 

lFMASf II » I=1,NI Jl 

GOGO 

58 

44 


CALL  PCHFCK 

IF^NZ,EQ,ICTR>  GO  TC-  IOC 


GO  GO 

C1LL  PSHEC<  TO  RESIN  PROCESSING  THE  P.RTICLE  OPT.  INTO  .  »»  gj 

SS,!S!U«,  so  t  c*  ,«o  SB 

rn  TAPE  FOR  USE  IN  SUBSEQUENT  GOGO 

CALL  POUMP  TO  OUMP  PARTICLE  DATA  ON  TO  TAPE  FO  U  GOGO 

MAP  corf  LOAOS  gogo 

CALL  PQHP  fJfO 

GOTO  100  rnr-o 

1,00  RETURN  G0C'° 

END 


c/oooooooooooooooonoooooooooonr^oooooooooooooooonooooooooo 
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SUBROUTINE  LINKS  (NUMTAP)  LINKS 

LINKS 

THIS  PROGRAM  INITIALIZES  AND  WRITES . HE AOINGS  FOR  THE  OUTPUT  LINKS 

PROCESSOR.  THEN  IT  CALLS  THE  FIRST  PART  OF  THE  PARTICLE  ACTIVITY  LINKS 
MODULE  (PAMil  TO  PRECOMPUTE  DATA  USED  RY  THE  SECOND  PART  OF  THE  LINKS 
PARTICLE  ACTIVITY  HCOULF  WHICH  WILL  BE  CALLEO  DURING  THE  LINKS 

EXECUTION  OF  LINK9.  LINKS 


MAOA  VEPSION 
H.G.NORMENT 


CUTMAP 
PELTAX 
OGX,  OGY 

FMASS(I) 

FP<I) 

FSUM 

FH 

GRUFF 

IC(J> 

10(1*1 .GT.O 
IC(1S) .GT.O 
ICTR 
IH 

IV 

INC 

IPOUT 

ISOUT 

ISIN 

IRROR 

ITAB 

JC(J) 

JC(1»=1 
JC (11 =2 
KTR(I) 

MARRAY 

HASCHN 

MXREO 

NE 

NIJ 

NHAP 

NOL 


NOR 


NORD 


LINKS 

-  DESIGNED  TO  OPERATE  WITH  THE  OTM  LINKS 

JUNE  25»l971  LINKS 

LINKS 

»»»*»«*»****»»,  GLOSSARY  **********#****<m>*****¥*»**lINK8 

LINKS 

CLT-OFF  THRESHOLD  FOR  HAP  ORDINATE  VALUES  LINKS 

MAXIMUM  WIDTH  OF  A  CORE-LOAD  HAP  LINKS 

HAP  GRID  POINT  SEPARATION  DISTANCES  IN  THE  LINKS 

X  AND  Y  DIRECTIONS  LINKS 

FALLOUT  HASS  FRACTION  IN  EACH  SIZE  CLASS  LINKS 

TOTAL  RADIOACTIVITY  IN  EACH  SIZE  CLASS  LINKS 

SLH  CF  ALL  MAP  POINT  ORDINATES  LINKS 

FISSION  YIELD  LINK8 

A  COMBINED  GROUND  ROUGHNESS  AND  RADIATION  METER  LINKS 
RESPONSE  FACTOR  LINK8 

RLN  CONTROL  VARIABLES  LINKS 

NO  MAPS  ARE  TO  BF  PRODUCED  LTNK8 

PRINT  CONTENTS  OF  TAPE  IPOUT  LINKS 

SEE  C-OGO  GLOSSARY  LINKS 

PRINTER  DESCRIPTION—  NUMBER  OF  CHARCTERS/INCH  LINKS 
ACROSS  A  PAGE  OF  PRINTEO  OUTPUT  LTNK8 

PRINTER  DESCRIPTION—  NUMBER  OF  CHARCTERS/INCH  LINKS 
DOWN  A  PAGE  OF  PRINTED  OUTPUT  LINKS 

NUMBER  OF  MAP  OROINATE  COLUMNS  THAT  CAN  PE  LTNKS 

ACCOMODATED  BY  THE  PRINTER  PAPER  LINKS 

TAPE  ON  WHICH  PARTICLE  PARAMETERS  ARE  WRITTFN  BY  LINKS 
TME  OTM  (OPH  BINARY  INPUT)  LINKS 

SYSTEM  OUTPUT  TAPE  NUMBER  LINK8 

SYSTEM  INPUT  TAPE  NUMBER  LINK8 

ERROR  STOP  TRACE  WORD  LINKS 

NUMBER  OF  PARTICLE  SI7E  CLASSES  LINKS 

MAP  SPECIFICATION  CONTROL  VARIABLES  LINKS 

2  LINE  E  FORMAT  LINKS 

2  LINE  F11.3  FORMAT  LINKS 

SEE  PCHECK  GLOSSARY  LINKS 

PARTICLE  DATA  ARRAYS  DIMENSION  LINKS 

MASS  CHAIN  NUMBER  FOR  A  NREQ=10  REQUEST  LINKS 

MAXIMUM  NUMBER  OF  PROCESSING  REQUEST  TYPES  LINKS 

ALLOWED  FOR  IN  THE  CODE  LINKS 

SEE  PCHECK  GLOSSARY  LINKS 

NIMBER  OF  PARTICLE  DESCRIPTIONS  IN  THE  CURRENT  LINKS 
PARTICLE  BLOCK  LINKS 

MAXIMUM  NUMBER  OF  HAP  POINTS  IN  A  MAP  CORE  LOADLINKS 
SMALLEST  X  INDEX  OF  A  MAP  POINT  TO  THE  RIGHT  OF  LINKS 
TEE  l EFT  BOUNDARY  OF  THE  CONTRIBUTION  ELLIPSE  LINKS 
OF  A  DEPOSIT  INCREMENT  LINKS 

LARGEST  X  INDEX  OF  A  MAP  POINT  TO  THE  LEFT  OF  LINKS 
TWE  RIGHT  BOUNDARY  OF  THE  CONTRIBUTION  ELLIPSE  LINKS 
OF  A  DEPOSIT  INCREMENT  LINKS 

ROUTING  PARAMETER  FOR  PARTICLE  CONTRIBUTIONS  LINKS 


1 

2 

3 

4 

5 

6 
7 
3 
Q 

10 

11 

12 

13 

1** 

15 

16 
17 
16 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 
36 

39 

40 

41 

42 

43 

44 

45 

46 

47 
46 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 


46 
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AT  MAP  POINTS  -  - 


1  - 


2  - 


7  - 


NOX 


NPFQ 

NRO 

NST 

NTASK 

NUMTAP(  ) 

NXMAP 


AXIS  IN  A 


NYMAP 


AXIS  IN  A 


NZ 


OMAP(J) 
OPIO(  ) 
PACT  II) 
PSIZE(I) 
OCUT 


T1,T2 

TEXIT 

TIME , rENTER 

X, Y,Z0UT,7, 
SYOT.ROUT 

XMAX,XMIN 

YMAX,YMIN 

XI,  X2 


WITM  ) 
WFMAS(I) 


ZOEP 


time  OF  arrival  (NREQ=11) 

SMALLEST  PARTICLE  SIZE  (NREQ=13) 

TIME  OF  CESSATION  CNRE0=12) 

LARGEST  PARTICLE  SIZE  CNREQ=14) 
STRAIGHTFORWARD  ADOITION  OF  THE 
GAUSSIAN  DISTRIBUTED  QUANTITY  TO  EACH 
MAP  POINT  (NP£Q=2-10,  15-16) 

4  -  COUNT  OF  DEPOSIT  INCREMENTS  (NRE0=1) 
NLMBER  OF  GRID  POINTS  ALLOWED  IN  X  DIRECTION 
IK  A  CORE-LOAD  HAP 
COMPUTATION  OPTION  CODE 
A  COUNTER  FOR  MAP  REQUESTS 
TALLY  OF  PARTICLE  DATA  BLOCKS 
A  TALLY  OF  MAP  SPECIFICATIONS 
TAPE  NUMBER  ARRAY 
NUMBER  OF  MAP  POINTS  ON  THE  X 
LCAO 

NLMBER  OF  MAP  POINTS  ON  THE  Y 
LCAO 

NLMBEP  OF  MAP  CORE  LOADS  REQUIRED  IN  ADDITION 
TFE  FIRST 

TFE  MAP  ORDINATE  ARRAY 
OLTPUT  PROCESSOR  IDENTIFICATION 
PARTICLE  SIZE  CLASS  UPPER  BOUNDARY  DIAMETER 
PARTICLE  SI7E  CLASS  CENTRAL  DIAMETERS 
CUT-OFF  THRESHOLD  FOR  AN  INDIVIDUAL  OEPOSIT 
INCREMENT  CONTRIBUTION 
REQUEST  TIME  ARGUMENTS 

TIME  RELATIVE  TO  SHOT  TIME  CORRESPONDING  TO 
TIME  RELATIVE  TO  SHOT  TIME  CORRESPONDING  TO 
SXOT,  PARTICLE  DESCRIPTION  PARAMETERS CALL  INDEXED) 
,PS , FNAS  SEE  GLOSSARY  IN  DTM  REPORT 

MAXIMUM  AND  MINIMUM  X  COORDINATES  OF  THE  MAP 
MAXIMUM  AND  MINIMUM  Y  COORDINATES  OF  THE  MAP 
X  AXIS  BOUNORY  COORDINATES  OF  THE  CURRENT  MAP 
CCRE  LOAD 

DTM  IDENTIFICATION 

TCTAl  MASS  OF  FALLOUT  IN  EACH  PARTIGLE  SIZE 
CLASS/  GRUFF 

ALTITUDE  OF  THE  PARTICLE  DEPOSITION  PLANE 


59 

60 


LINKS 
LINKS 
LINKS  6). 
LINKS  6? 
LINK8  63 
LINKS  64 
LINK8  65 
LINKS 
LINKS 
LINKS 
LINKS 
LINKS 
LINKS 
LINKS 
LINKS 


66 

67 

68 
59 

70 

71 

72 

73 


LINKS  74 
MAP  CO RELINKS  75 
LINKS  76 
MAP  CO  RELINKS  77 
LINKS  78 
T0LINK8  79 
LINKS 


GLOSSARY  OF  CCMMCN  /SETl/  SEE  OASA-1SOO-III (REVISED) 


MMON  7SET17 

»OETIO(12 )  ,OIAM(201) 
ID (200) , IDISTR  , IEXEC 
;tr  .TTnt?nn» 


J  AV  4.  — »  « 

»TIO (200 ) 
,  T2M 

3HB  , NHODO 

ION  /PAROAT/ 

10)  ,Y(500) 

,FMAS(500 


, IEXEC 
,S0 
»U 

,7V(200) 


,OMEAN 
» IRTSE 
»  SSAM 
,  VPR 

, VX (200) 


»ONS 
,ISIN 
,TME 
♦  H 

,VY (200) 


,EXPO 

,ISO'JT 

»TMP1 

,HBURST 


,ZOUT (500) 
»  ,KTR(500I 


,SXOT (500) 
»  F 


,SYOT (500)  ,ROUT (500) 
,  GAMA  ,  BSO 


S3 

84 


T2 

T1 


37 

88 

89 


SO 

LINK8  81 
LINKS  82 
LINKS 
LINK8 
LINKS  85 
LINKS  86 
LINKS 
LINKS 
LINKS 
LINKS  90 
LINKS  91 
LINKS  92 
LINKS 
LINKS 
LINKS  95 
LINKS  96 
LINKS  97 
LINKS  98 
LINKS  99 
LINK8100 


93 

94 


►LINKS101 


LINK8102 
►♦LINK8103 
LINK8104 
LINK8105 
LTNK8106 
L1NK8107 
LINKS10S 
, LINK8109 
, LINKS110 
, LINK8111 
,LINK8112 
LINK8113 
LINK6114 
, LINK8115 
, LINK6116 
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■M 


1 


-7,^- 

- 

3ASQ  ,SINA 

,  COSA 

,HFMAS(200) 

,YPRMU 

, YPRML 

jLtNKSllZ 

4TC500) 

LINK0118 

COMMON  /RUNOAT/ 

LTNK8119 

INI J  ,NE 

,NPEQ 

}N7 

*ICTR 

,NXMAP 

, LINK8120 

2T1  ,T2 

,,MAPr<UN 

*  TG2 

tl° 

, 1C (19) 

,LIN<8121 

3ICC18)  ,NYM«P 

,NTAS< 

,  NORD 

,XGZ 

»YGZ 

LINK8122 

COMMON  /MAPOAT/ 

LINY3123 

10M AP f 15000 ) ,QCUT 

,CUTMAP 

»OGX 

♦  OGY 

,DELTAX 

, LTNK8124 

2XMAX  , XMIN 

,  YMAX 

,  Y  MIN 

,FSUM 

,RUFSAH 

,LTNK8125 

3X1  ,X2 

, MFTAPE 

*  ZOEP 

LINK8126 

COMMON  /CONDAT/ 

LIHK8127 

1I°0'JT  ,  JPOUT 

,K.  .  JT 

,  KT  APE 

,LT  APE 

,MARRAY 

,LTN<8128 

2NM A°  , MXREO 

»IH 

,IV 

LINKS 12 9 

COMMON/OUTPUT/ 

LINK3130 

1  FISNUM  , 

FP  (200) 

*  FW  , 

ITAB 

» JGO 

LINK8131 

2  » MASCHN  , 

PS12E  <20  0) 

,FMA5S 1200)  , 

PAC-TI200I 

LINK813? 

LTN<6133 

C 

*** *******  ********************  ***<t****t.  *********** *■**•*• *************  LI NK8134 

r* 

V 

LINK8135 

c 

DIMENSIONS  PECULIAR  TO  LINK* 

LINK8136 

DIMENSION  CRIO (12) 

♦0PIDC12)  ,PSEin<12i 

LINK8137 

1 , NUMT  AP (15 ) ,WIO '12) 

LINK8138 

OlMEMStON  00(260 

LTNK8139 

C 

LINK8140 

c 

►’••♦LINK8141 

c 

LINK8142 

1 

FORMAT (12 A6) 

LINK8143 

2 

FOPMAT (15X , 181 4) 

LINK9144 

to  F0PMAT(////29X,63H****  SUMMARY  OF  PRECEDING  DELFIC  MODULE  RUN  IGELIN<«145 
1NTIFIE9S  ****  III 25Xh3H***»  OUTPUT  PROCESSOR  IDENTIFICATION  *»*LTNK8146 
2*/25X,12A6///25X,56H'«**  INITIAL  CONDITIONS  I  FIREBALL)  IDENTIFICALIN<814? 
3TION  ****/25X, lCu6///25X,37rt***»  CLOUD  RISE  IDENTIFICATION  ****LINK8148 
4/20 ,12A6///25X,57H****  CLOUD  RISE-TRANSPORT  INTERFACE  IOENT JFICALINK8149 
5TION  *•**,  725X  , l2A6/7/25X»46H*  ***  DIFFUSIVE  TRANSPORT  IOFNT IFLINK8I50 


6ICATICN  ***-/25X, 12A61  LINK9151 

11  FORM,*T(//15X2t  HTRANSPORt  IOENTIFIC  ATION/'725X»  12A6-I  LINK8152 

12  FORMAT J//25X.24H****  OTHER  INPUTS  *•**)  LINK8153 

15  FORMAT (1814)  LIN<8154 

16  FORMAT (/15X77H****  THE  CONTROL  VARIABLE  ARRAY,  IC( J) ,  HAS  GIVEN  TLINK8155 

1HE  FOLLOHING  VALUES  *»**)  LINK8156 

21  FORMAT (//15X43HPRINTER  DESCRIPTION  -  CHARACTERS  PER  INCH)  LINK8157 

22  FORMAT  tlSX , 14HHORIZCNT  ALI5 , 1 OX, 3  CU VERTICAL  III  uTN<8158 

26  FORMAT C17X,1HX,11X,1HY,  llXj-'HZ*  11X,1HT*  9X,4HSX0T9  8X,4HSY0T,  LINK6159 

1  8X,  4HR0UT,  9X,2HPS,  9/-4HFHmS/M  LINK8160 

28  FORMAT (1H1///51X19H*  ********  *//12X101HT  HE  D  E  P  A  R  TLINK8161 
1MEST  OF  DEFENSE  FALLOUT  PREOICTI  0LINK8162 

2  H  S  Y  S  f  E  M,//51X,19H*  ********  *////48X , 23HOUTPUT  PR0LINK8163 

3CESSOR  M00UL"///55X,11HPREPARE0  BY744X,31HHT.  AUBURN  RESEARCH  ASS0LINK8164 
4C.  ,INC../54X,13HNEHTCN,  MASS.)  11(^:155 

29  FORMAT t;////45X41HLISTING  OF  DEPOSIT  INCREMENT  DESCRIPTIONS)  LINK8166 

3G  eORMAT (7/13X6HBL0CK  :4)  LINK3167 

36  FCPMSTC10X,9E12.4)  LINK8168 

37  FORMAT  (lQIt ,  43HN0.  CF  DEPOSIT  INCREMENTS  IN  THIS  BLOCK  IS  I4i  LINK8169 

39  FORMAT (46H  NO  MAPS.  THIS  RUN  FOR  TAPE  IPOUT  PRINT  ONLY.)  LINK8170 

C  LINK8171 

r*********************,c**  BEGINNING  OF  PROGRAM  *************** ******* **llNKilf 2 

C  LINK8173 

LOGICAL  SKIP  LINK8174 
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data  progrk  /sh  links/ 
kout=isout 
IP0UT=NUMTAP(9> 
jpnor=NUMrap(  2* 

KP0'JT  =  NUMTAP(  3) 

MXPEP=20 
marray  =  500 
NMAP=1<5C00 
NTASK=0 

reao  ipour  header  data 

REWINO  I POUT 
PE  AO  (IPOUTUPOTJ 

READ (IPOUT I FW,SSAM,SL0TMP,TMSD,SIGHA,TW,HBURST,X6Z,YGZ,TGZ,9Z, 
lROPaRT,IRAn,RADM&X,ZBRSTZ 
RT  MRURST  in  METERS  TO  HOP  IN  FEET 
HOP=H0'JPST  A.  3  0  48 

REao  PREVIOUS  IDENTIFIERS  FROM  GPOUNOEO  PARTICLES  TAPE 
READ ( IPOUT ) (OSEIDT J),J=1,12> 

RE  AO ( IPOUT) I  CRIO( J), J=l,12) 

REaodPOUTI  (OETIO(JI,J=l,12> 

RFAO { I cOUTT (WIO(JT , J=1.121 
READ  (IP0UT1ITA0 

REaOTIPOUTI (PSIZE(J),FMASS(J) ,PACT(J» ,J=1»ITA9I 
REAO (IPOUT • NAT 

REAOCIPCun  nD(Jl,00(J>,D9<iJi  ,  J-l,  NAT) 

READ  IDENTIFIER  FOR  OUTPUT  PROCESSOR  RUN 
REAO  (ISIN,1> (OPIO(  J)  ,  J=l,12» 

REAO  CONTROL  VARIABLE  ARRAY 
PE  AO  (ISIN,15»  <IC( J), J=l,18l 

THIS  PART  OF  THE  COTE  DUMPS  TAPE  IPOUT  IF  REQUIRED 
IC  < 18 1  POSITIVE  MEANS  OUMP  TAPE  IPCUT  BEFORE  EXECUTION 
SKIP  =  .TRUE. 

IC(lR)  =  0  MEANS  DO  NOT  DUN®  TAPE  IPOUT 
IF( IC (181 I  500,5021,502 
500  IRROR=-500 
GO  TO  333 

502  SKIP-. FALSE. 

WRITE  IISOUT, 281 

WRITE  ( ISOUT ,11)  C'lCTJ* ,J=1,12» 

WRITE  (ISOUT, 291 
5021  NST  =  0 
600  READ  (IPOUTVNIJ 
NST=NST*1 

IF(NI%I)  503  ,50  1,50  A 

503  IOR0R=-503 
GO  TO  333 

504  READ  (IPOUT  I  (X(I),T  (I)  ,ZCUT<I>  ,1  (I)  ,SXOT(I»,Srr,m»,ROUT(I),PSi  II  , 
lFMAS(I) ,I=1,NIJI 
IF (SKIP)  GO  TO  600 
WRITE  (ISOUT, 30INST 
WRITE  (ISOUT, 3  Tl NI J 
WRITE  (ISOUT, 261 

WRITE (ISOUT, 361 ( X < I I ,V ( Ti , *OUT ( II . T ( I > ,SXCT(I» ,SfOT (II ,ROUT (I) , 
IPS  (I)  .FMASIII  ,I=1,NIJ) 
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LINK3176 
LI  K‘KB177 
LINK517  9 
LINKA179 
LTNK3130 
LINKB1M 
LTNK313? 
LINK81B3 
LINK3134 
LINK8185 
LTNK8186 
LINK5187 
LINK8188 
LTNK31S9 
LINK8190 
LINK3191 
LINK8192 
LINK8193 
LINK8194 
Li NK8195 
LI  '1K319F 
LINK8197 
LINK8196 
LINK8199 
LINK8200 
LINK8201 
LTN’K8202 
LTNK8203 
LINK8204 
LINK3205 
LINK8206 
LINK3207 
LINK8208 
LTNK8209 
LINK8210 
LINK8211 
LTNK8212 
LINK6213 
LINK3214 
LINK8215 
LINK321F. 
LINK8217 
LINK8218 
I.INK8219 
LINKS220 
LTNKB221 
LINKR222 
LINK3223 
LINKR224 
LIN<8225 
LINK8226 
LINK8227 
LINK822C 
LINK8229 
LINK8230 
LINK8231 
LINK8232 
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Ul)  O  O 


501 


505 

506 

333 

510 


GO  TO  600 
CONTINUE 

IC( 17)  POSITIVE  MEANS  STOP  WITHOUT  ENTERING  OUTPUT  PROCESSOR 
IC(17)  =  0  WEANS  PROCEED  WITH  JOft 
IF  C IC (17) I  506,511,510 
IRPOR=-506 

CALL  ERROR  (PRGGRM, IRROR»ISOUT) 

WRITE  (7SOUT,  39) 

CALL  EXIT 

ENO  OF  TAPE  IPQUT  OIMP 


J)»J=1»12)»(CRID(J>»J= 


II  CONTINUE 

PE AO  PRINTER  DESCRIPTION  -  CHAR/INCH  HORIZON! A L, VERTICAL 

III  REAO  ( ISIN ,15) IH » I V 

PRINT  A  HEADING  TO  IOENTIFY  PRINTED  OUTPUT 
WRITE  (IS0UT,28) 

WRITE  ( ISOUT, 10)  (OFID(J) » J=1 » 12) * (DET ID( 

12), (PSEIO( J),J=1,12>, (WIO( J) , J=l,12» 

WRITE  (ISOUT, 121 
WRITE  (ISOUT* 16) 

'.'RITE  (ISOUT, 2)  CIC(J),J=1,18) 

WRITE  (ISOUT, 21) 

WRITE  (ISOUT, 22) IH, IV 

CALL  PAM1 

1  (HCB  , SLDTHP  , TMSO  ,TW 

2  , ISIN  , ISOUT  , IPCUT  ,NUNTAP  , SIGMA  I 
RETURN 
ENO 


LIN<*233 
LINK6234 
LINKA235 
LINK623C 
LTN<8237 
LINK823R 
LINK9239 
LlN<8240 
LINKB24) 
LTN<8242 
LINK8243 
LINK9244 
LIN<«245 
LINK8246 
LlN<824  7 
LTNW8248 
LTN<8249 
1 ,1LINK8250 
LIN<8251 
LINK825? 
LINK325? 
LIN<8254 
LINK8255 
LIN<8256 
LINK8257 
LINK3259 
LIN<8259 
LINK8260 
LINXL8261 
LINK8262 
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SUBROUTINE  LINK9 

MAOft  VERSION  -  OESIGNEO  TO  OPERATE  WITH  THE  OTM 
H.G.NORHENT  JUNE  28,1971 


LINK9  1 
LTN<9  2 
LINK9  3 
LINK9  4 
LIN<9  5 

LlNiCR  7 

SECONO  HALE  OF  THE  CUTPUT  PROCESSOR  MAIN  CONTROL  PROGRAH  LINK9  6 

THIS  SUBROUTINE  INITIALIZES  AND  CONTROLS  FOR  HAP  CALCULATIONS  LTN<9  9 

LINK9  10 
•LINK9  11 
LINrfR  12 
LINK9  13 
LlN<9  14 
®*LINK9  15 
LTN<9  16 
LINK9  17 
LINK9  lft 
**LIN<9  19 


SUBROUTINES  CALLED  - 

GOGO 

HAP 


GLOSSARY  •**»*♦*•***♦***»*♦♦*****» 

SEE  THE  LlN<8  GLOSSARY 


LTNK9 

20 

COMMON  /SET  17 

LlN<9 

21 

1CAY  , CE  T ID ( 1 2  ) 

,niAM(201) 

.DMEAN 

,ONS 

♦  EXPO 

• LINK9 

22 

2OITIO{2O01 , IOISTR 

,IFXEC 

, IRTSE 

,  IS  IN 

,ISOUT 

, LIN<9 

23 

3NOSTR  , TIO (200) 

,  SO 

,ssam 

,TME 

.tmpi 

, LINK9 

24 

4TMP2  ,T2M 

»U 

,  VRR 

»w 

,H9URST 

, LINK9 

25 

5SCL0HB  , NHOOC 

,ZV(20Q) 

,VX(200) 

,VY«200) 

LINK9 

26 

COMMON  /PAROAT/ 

LINK9 

27 

1X(F(50)  ,  Y(500) 

,ZOUT (500 ) 

,SXOT (500) 

,SYOT (500) 

,ROUT (500) 

,LIN<9 

28 

2PS (5001  ,  FMAS  <50  0  ) 

,<TR(500) 

»  F 

,  GAMA 

,  BSQ 

, LINK9 

29 

3ASO  «  SIN A 

,COSA 

,WFMAS(?00) 

,YPRMU 

,  YRRML 

, LIN<9 

30 

41(500) 

LINK9 

31 

CHMMON  /PUNOAT/ 

LIN<9 

32 

1NIJ  , NE 

,NREQ 

,N2 

,IC  TR 

,NXMAP 

, LINK9 

33 

2T1  ,T2 

, MAPRUN 

,  T  GZ 

,IP 

,  JC(13) 

, LINK9 

34 

TIC (18)  ,NYMAP 

,NTAS< 

,NORO 

,XGZ 

♦  YGZ 

LIN<9 

35 

COMMON  /HAPGAT/ 

LINK9 

36 

1QMAPC15000* ,QCUT 

,CUTMAP 

,OGX 

,OGY 

tOELTAX 

, LIN<9 

37 

2XMAX  ,XMIN 

,  YMAX 

,  YMIN 

,FSUM 

♦RUFSAM 

, LINK9 

32 

3X1  ,X2 

, MPTAPE 

,ZDEP 

LINK9 

39 

COMMON  /CONOAT/ 

LINK9 

40 

1IPOUT  , JPOUT 

, K°OUT 

,KTAPE 

,LT  APE 

♦MARRAY 

, LINK9 

41 

2NMAP  , MX RED 

»IH 

,IV 

LlN<9 

42 

COMHCN/OUTPUT/ 

LINK9 

43 

1  FISNUM  , FP 

1200)  ,FU 

1  ? 

IT  A8 

♦  JGO 

LIN<9 

44 

2  , MASC.HN  ,PSIZE  T20 O'  ,FMASS(200)  , 

PACTI200) 

LINK9 

45 

COMMON/OECAY/ 

LIN<9 

46 

1  IGO  ,J0 

f <OOS  f 

TENTER 

LINK9 

47 

2  , TEXIT  ,TIME 

LINK9 

48 

LINK9 

49 

••LINK9 

50 

LIN*9 

51 

F0PMAT(12A6> 

LINK9 

52 

FORMAT  (//15X, 23HSUM 

OF  MAP  ORDINATES  =  E13 

.6  ) 

LIN<9 

53 

FORMAT (1H1////54X , 1 1H*  *  *  *  *  » 

) 

LINK9 

54 

FORMAT  (//15X,  23HGRCu::0  ROUGHNESS  FACTOR  F10. 3,10X,15HALTITUOE  OF 
1GZ  F10.3) 

FORMAT  (7F10.3) 

FORMAT  (18T4I 


L:Nk9  55 
LIN<9  56 
LINK9  57 
LlN<9  58 
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17  FOPMAT  (32H  OUTPUT  PROCESSING  IS  COMPLETED.! 

23  FORMAT (1H1///39X27H****  OUTPUT  PROCESSOR  TAi:<I5,6H  **•*> 

24  FORMAT C///15X25HGRI0  LIMITS  AND  INTERVALS/18X4HXHIN10X4HXHAX1QX4 
1HIN1GX4HYMAX10  X7H0ELT A  X1GX7HDELTA  Y/15XF10.0 ,4XF10.0 ,4X^10. 0 ,4X 
20.0,5XF10.1t5XF10.1T 

29  FORMAT I//15X71HTHE  CONTROL  VARIABLE  ARRAY,  JC(J),  HAS  BEEN  GIVEN 
1HE  FOLLOWING  VALUES ./15X18I4) 

27  FORMAT (715X32HMAPPEC  ON  GRID  INTERVALS  DGX  =  P10.i,7H  OGY=FlO 
1) 

31  FORMAT (117HQlNAO£GiU  ATr  PRINTER  DESCRIPTION.  AN  UNOISTORTEO  NAP  C 
1NOT  8E  GUARANTEED.  THIS  RUN  WAS  CONTINUED  WITH  PRINTER  DESCRP./5 
214WIH=1 0  ANC  I  V=6) 

32  FORMAT(I4,2F10.3,I4,2F10.3> 

33  FORMAT (25H0UNACCEPTA8LE  REOUEST  ...14) 

34  F0RMAT(//////15X,15FRECUEST  NUMQER  14/ //15X,5HTYPE  I4,10X5HTl  ^ 
10. 1 »1QX»5HT2  =  F10. 1,19X,9HMASCHN  =  14/7  15X,6HQCUT=  ,E12.5,10X, 
2CUTHAP=  ,E12«  51 

35  FORMAT ( 15) 


DIMENSION  nn<?60) 

OATA  PROGRM  /&H  LINK97 
LOGICAL  JD,<O0S,IGO 

IGO=.TRUE. 

NUL=0 
FSUM=0.0 

119  IF  (  FSUM) 1601, 1191, 1601 

1601  WRITE  f ISOUT ,21  FSUM 
FSUM-0.0 

1191  READ  LIMITS  ON  AREA  CF  INTEREST 
1191  READ  (ISIN, 9IXMAX , XH IN, YMAX , YMIN,OGX,DGY , GRUFF 
IFC GRUFF) 1602, 1602, 1603 

1602  GRUFF=1.0 

1603  IF (A8S (OGX)  «•  A8S(0GY1 1120,120,121 
120  WRITE  (ISOUT, 17) 

REWIND  I POUT 
PRINT  17 
RETURN 

READ  OTHER  SPECIFIC  INPUT 
121  REA0IISIN,9)Z0EP 

READ  <ISIN,15) (JC( J),J=1,1») 

NTASK=NTASK*1 
NRQ=0 

CHECX  PRINTER  DESCRIPTIONS 
NI=IH*IV 

IF(NI)601,601,122 
601  IH-10 
IV=6 

WRITE  (ISOUT, 31) 

C 

C  122  WRITE  A  LOCAL  HEACING 
1 22  WRITE  (ISOUT, 23)NTAS< 


LINK9  59 
LIM<9  60 
HYLTNK9  61 
F1LINX9  62 
LINK9  63 
TLTN<9  64 
LINK9  65 
. 1LINK9  66 
LIN<9  67 
ANLINX9  68 
X,LIN<9  6« 
LINK9  70 
LTNK9  71 
LINKO  72 
F1LINK9  73 
8HLINK9  74 
LINKO  76 
LIN<9  76 
LINK9  77 
**LTN<9  76 
**LINK9  79 
LTN<9  60 
LINK9  81 
LINK9  92 
LtN<9  83 
LINK9  84 
LTN<9  85 
LINK9  86 
LINK9  87 
LINK9  88 
LINK9  99 
LIN<9  90 
LINK9  91 
LINK9  92 
LIMK9  93 
LTNX9  94 
LINK9  95 
LINK9  96 
LINK9  97 
LINK9  98 
LTN<9  99 
LINK9100 
LINK9101 
LINX9102 
LINX9103 
LINX9104 
LINX9105 
LTNK9106 
LINK9107 
LINK91Q9 
LIN<9109 
LINK9110 
LINX9111 
LTNK911? 
LtN<9113 
LINK7114 
LINK9115 
LINX9116 
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1225  HP ITE  ( ISOUT ,24  I  XMIN,XHAX, YHIN,YMA X,OGX,DGY 
W°ITE  (ISOUT, 4)  GRLFF  ,ZOEP 
WRITE  (IS0UT,25> CJCCJI ,J=1,H) 

GO  TO  1209 


C 

1211 

C 

1209 
1604 
C  121 
1219 


CONTINUE 


TF(FSUM) 1604,  1219, 1604 
WRITF  ( ISOUT, 2)  FSUK 
9  READ  A  REOIJEST  FOR  PROCESSING 
REAn(ISIN,32»NPEQ,Tl,T?,MASCHN,0CUT,CUTMAP 
IF(MASCKN.EC.0.ANO.»'REC-.NE.10>GO  TO  1210 
:P(HASCHN.GT.ri. ANO.NASCMN.lt. 1621  GO  TO  1210 
CALL  ERRORCPROGRrt, 1209, ISOUT) 

MA$CHN=95 


C 

C 

C 


CLEAR  OUT  THE  OHAP  ARRAY 


1210 


935 


CLROT=0.0 

IF{ (NREQ.EQ.il) .OR.  1NRE0.E3. 13) )  CLROT=l.E30 

no  035  I=1,N<AP 

OHAP  f I ) =CLROT 

M A  PRUN  =  0 

NRQ=NRQ4-1 

FSUM  ^  0.0 


C 

C 

C 


1213 

1215 


IS  NREQ  AN  ACCEPT A9LE  REQUEST 
NO  TO  1215 

Ic(NREO) 1212,1212, 1213 
IF {NPEO-MXPEO) 1214,1214,1215 
WRITE  ( ISOUT, 33) NREC 
GO  TO  1209 


C 

C 

C  121 
1212 


MXPFO  IS  MAXIMUM  NUMBER  OF  CALCULATION  COOES  ALLOWED  FOR  IN  CALC. 
2  NO  MORE  REQUESTS.  PREPARE  TO  RETURN  TO  READ  LOCAL  DATA. 

CONTINUF 
GO  TO  119 
1214  REWINO  TPQUT 

READ(IPOUT)  ITST 

READCIPCUT)TST,TSTt  TS T , TST , TST, TS7 , TST, TST, TST , TST, TST, TS T, I TST, 

1  1ST, TST 

READ(IPOUT) tOETIOT J), J=l,12) 

READ  ( IPOUT) (DETIO( J), J=l,12) 

RE  AD  * IPOUT ) COETIOl J), J  =  l,12) 

READ(IPOUT)  (OETID()> ,J  =  1,12J 
READ (IPOUT) IT  AP 

READ (IPOUT) (OD(J)»OC(J),ODf J),J=1,  I  TAB) 

READ (IPOUT) NAT 

READ! IPOUT) (ODIJ) ,OC(J) ,OD( J) , J=1,NAT) 

IF (NRO-1) 1221, 1222, 1221 
WRITE  (ISOUT -  3) 

WRITE  (IS0UT,34)NRQ,NPE0,T1,T2,MASCWN  ,OCUT,CUTMAP 
RUFSAM  =  SSAH7GRUFF 
JGO=  l 

FISNUM=FW*1.45E15 

IF(NREQ.EQ.15.09.NREO.E0.16)  GO  TO  1223 
T1=TIMSEC(T1,3,0) 

T2=TIMSEC(T2,3,0) 


1221 

1222 
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LINK9124 
LIN<9125 
LTNK912F 
LIN<912  7 
LIN*9126 
LTN<9129 
LINK913  0 
LINK9131 
LIN<913? 
LINK9133 
LTN<9134 
LINK9135 
LTN<913E 
LINK9137 
LINK9135 
LIN<913Q 
LTNK9140 
LIN<9141 
LINK9142 
LIN<9143 
LINK9144 
LINK9145 
LINK9146 
LINK9147 
LINK914P 
LINK9149 
LIN  <9150 
LINK9151 
LINK9152 
LINK9153 
LINK9154 
LIN<9155 
LINK9156 
LTN<9157 
LINK9156 
LTNK9159 
LlN<qi60 
LINK9161 
LINK9162 
LINK9163 
LTN<9164 
LINK9165 
LTNK9166 
LINK9167 
LINK9168 
LIN<9169 
LINK9170 
LINK9171 
LINK9172 
LINK9173 
LINK9174 
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1223  TIME=Tl-TGZ 


r£NTER=TIHE 

LINK9175 

TEXIT=T2-TGZ 

1INK9176 

JO=. TRUE, 

LINK9177 

<OOS=, FALSE, 

LINK9179 

31 

LIN<917Q 

10LINK9190 

TIME=3600»  0 

LIN<9181 

32 

GO  TO  7Q 

LINK9182 

<0OS=. TRUE  « 

IIN<9163 

13 

GO  TO  78 

LINK9184 

JG0=2 

LTNK9195 

FXSNUP=FISNUM*l.Ef4 

LIN<9196 

rufsah=ssam 

LTNX919 7 

78 

GO  TO  Y9 

LIN<9188 

jo=. false. 

LINK9199 

79 

FISNL'M=FlSNUH/3600. 

LINX9190 

CALL  PAM2 

LIN*9191 

80 

continue 

LINK919? 

UNX9193 

INITIALIZE  FOR  PROCESSING 

LINK9194 

1705 

00  1705  I=1,ITAB 

LIN<9195 

MFHAfCII  -  FMASSCI) *RUFSAM 

LINK9196 

801 

IF  (NR0.GT.1)  GO  TO  1502 

LtN<9197 

LINX9198 

1217 

CONTINUE 

LINX9199 

1302 

dtsy=iv 

LINKR20C 

Rn=2.o*oiSY/orsx 

DTSX=RO*DGY 

.  _  ifrfnisx-DGxii3o3,i3C4,i304 

1303  0GX=0ISX 

1304  DGY=OGX/RD 

c 

c  PREPARE  TO  PROCESS  CUTPUT 

i«  ^H«4pc=L;^»r-i«",^ omc  toaos  re™™  f^RST 

NXMAP= ( XMA  X-XMIN ) /OCX 
NST  =  NMAP/NYMAP 

IF<NKHAP-NSn  1503,1503,1402 
1402  NXMAP=NST 
XMAP=NXM«o 

Z7=|XHAX-XMIN) /IXMAF*DGX) 

NZ=ZZ 

TST=NZ 

<en„  *F(ZZ'Tsn  1500,1501,1401 

1500  IRROR=-i5Q0 
GO  VO  333 

1501  N7=NZ-1 

GO  TC  1401 
1503  NZ=0 

'IZ'ZJSR'ST  "P  C0RE  Loaos  CALCUUri°” 

IF CN0X11403, 1403,14 C4 
1403  IRR0R=-14Q3 

CALI.  ERROR  (OR0GRH  ,  IROOP,  ISOUH 
GO  TO  1211 


333 
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LINK9201 

l!N<920? 

LINK9203 

LlN<9204 

LINK9205 

LIN<9206 

LINK9207 

LTNK920O 

LINX9209 

tIN<9210 

UNK92H 

LINK9212 

LINK9213 

LTNK9214 

LIN<9215 

LINK9216 

LINK9217 

LINK9215 

LIN<9219 

LINK9220 

LIN<9221 

LINK9222 

LTNK9223 

LINK9224 

LINK9225 

LINX9226 

LINK9227 

LTN<9228 

LINK9229 

LINK9230 

LIHK9231 

LINK9232 


haVi-rar  '&‘±£r^r. 


«@rsHs*«i*isl 


1404  OX=NOX 

dfltax=ox*ogx 

1502  WRITE  <ISOUT,27»OGX,OGY 

X1=XMIN 

X?=X1+0ELTAX 


300  CONTINUE 

301  TCTR=0 

IF(NZ) 203,204,207 
203  ie?or)0=-203 

GO  TO  333 

;  thts  is  the  branch  fop  a  single  core  load  map 
% 

✓ 

204  ktapf=ipout 
CALL  GOGO 
PEWINO  KTAPE 

1F(  (NREO.NE.llJ  .ANO.;.4REO.NE.13»)  GO  TO  305 
OO  302  IMAP=1,NMAP 

TF(OMAP(IMAP».G£.1.E30>  OMAP(IMAP) =0.0 

302  CONTINUE 

305  CALL  MAP 
GO  TO  1211 

3  THIS  IS  THE  BRANCH  FOP  A  MULTIPLE  CORE  LOAD  MAP 

207  <TAPF=IPOUT 
LTAPE= JPOIJT 
CALL  GC’GO 
REWIND  KTAPE 
HRITE(LT  APE1NUL 
REWIND  ltape 

IF(<NREQ.ME.lll.ANO.tNBFO.NE.13!>  GO  TO  308 
00  306  IMAP=1,NMAP 

IFTQMAOTIMAPT .GE.1.E30)  OMAPTIMAPI =0.0 

306  CONTINUE 
308  CALL  MAP 

00  220  IN0EX=1,N7 

C 

;  clear  out  tj€  omap  array 

c 

CLROT-O.O 

IF ( (NREQ.EQ.il) .OR. TNREO.EQ.  1 3)  I  CLROT=1.E30 
00  702  IHAP=1,NMAP 
702  OMAP(I)=CLROT 

IF(MC0(IN0EX,2).E0.1)  C-0  TO  208 
K1 APE— KPOUT 

ltape=jpout 
GO  TO  209 

208  <T APE= JPOUT 
LTAPE=KPOUT 

209  ICTR=INOEX 
X1  =  X2 

X2=Xl*0FLTAx 


LINK923? 
LINK9234 
LINX9235 
LTNK9236 
LINX9237 
LINK9238 
LINK9239 
LINK9240 
* LTNK9241 
LINX9242 
LTNK9243 
LTN<92'»4 
LINK9245 
LINK924F 
LINX9247 
LTNK9248 
LlN<924Q 
LINK9250 
LTNK9251 
LINK9252 
LTNK9253 
LIN<9254 
LINK9255 
LINKP256 
LINK9257 
LINK9258 
LINK9259 
LTN<9260 
LINK9261 
LINK9262 
LINK9263 
LINK9264 
LIN<9265 
LINK9266 
LINX9267 
LINK9268 
LINK9269 
LINK9270 
LINK9271 
LIN<9272 
LTNK9273 
LINK9274 
LINK9275 
LINK9276 
LIN<9277 
LINK9278 
LINK9779 
LINK9280 
LINK9281 
LINK9282 
LINK9283 
LINK9284 
LINK9295 
LINK9286 
LTNK9287 
LINK9288 
LINK9289 
LINKP290 


CALL  GOGO  LINK9291 

REHINO  KTAPE  LIN<9292 

WRITE (LT APE INUL  LINK9293 

REWINO  LTAPE  LINK929A 

IF ( ( NREO»NE»ll )»AN0.INRFQ.NE.13))  GO  TO  220  LIN<9?95 

00  215  IMAP=1,NMAP  LINK9?9f 

IF«OMAP(IMAP).GE.1.E30)  OMAP l IH&Pl =0. 0  LIN<9?:97 

215  CONTINUE  tTNK9296 

220  CALL  MAP  IIN<9?99 

GO  TO  1211  LINK930G 

END  LTNK9301 
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***MAP 

q 

MAF 

10 

COMMON  /SET  1/ 

MAP 

11 

1CAY  ,OEriO(12»  ,riAM<201)  , DMEAN 

,ONS 

}  EXPO 

,MAP 

12 

2DTTIO<2P0l  »IDISTR  , IEXEC  ,IRISE 

,  ISIN 

, ISOUT 

,  MAP 

13 

3NOSTR  , TID  (200 )  ,St>  ,SSAM 

,  THE 

,  TMP1 

,  MA  P 

14 

4TMP2  , T2M  ,U  ,VPR 

tM 

,HBURST 

,  MAP 

15 

5SCL0HR  , NHOOO  ,ZV<200)  ,VX(20Q» 

,VY  (200) 

MAP 

16 

COMMON  /RUNOAT/ 

MAP 

17 

INI J  ,NE  ,NREO  ,NZ 

,ICTR 

, NXMAP 

,MAP 

18 

2T1  ,T2  , MAPRUN  ,TGZ 

»IP 

,  JC  (1 8  ) 

,  MAP 

19 

3IC (18)  , NY  M  A  P  ,NTAS<  ,NORD 

,  XGZ 

,YGZ 

MAP 

20 

COMMON  /MAPOAT/ 

MAP 

21 

10MAR(19000» ,QCUT  ,CUTMAP  ,DGX 

,DG  Y 

,DELTAX 

,  MAP 

2? 

2XMAX  ,  XMIN  , YMAX  , YMIN 

,FSUM 

,rufsam 

,  MAP 

23 

3X1  ,X2  , MPTAPE  , 70EP 

MAP 

24 

COMMON/OUTPUT/ 

MAP 

25 

1  FISNUM  ■)  FP  <20  0)  ,FW 

,  ITAB 

,  JGO 

MAP 

26 

2  , MASCHN  , PS IZE  <200)  ,FMASS(200) 

,  PACT  <  2  00) 

MAP 

27 

MAP 

28 

♦  » 

*  **MAP 

29 

MAP 

30 

DIMENSION  JMAP (20 ) 

MAP 

31 

INTEGER  3LAN< 

MAP 

32 

DIMENSION  FMTEXP(21),FMTPUT(21) 

MAP 

33 

DATA  FMTEXP(l) , FMTRtT (1),FMTEXP<21), FMTRUTC21 ) , BLANK, 

FMT  A ,FMTF , 

MAP 

34 

1  FMTI/6H  (/IX ,  ,6H5X,  ,6H»  ,6H) 

,6H 

, 6HA  6  , 

MAP 

35 

2  6HR6.3  ,6HI6  /.00T/6H  .  f 

MAP 

36 

MAP 

37 

DATA  BITLUM,lNCfLRE)</  6HMULTIB,19,  0/ 

MAP 

38 

MAP 

39 

** 

***HAP 

40 

MAP 

41 

1 

FORMAT <1H1,5HSTRIPI 2) 

MAP 

42 

2 

FORMAT (/IX, 1916) 

MAP 

43 

3 

FORMAT ( 15X21HTW0-L INE  E  FORMAT  MAP) 

MAP 

44 

4 

FORMAT (5X,19F6. 3) 

MAP 

45 

5 

FORMAT (15X26HTW0-LINE  F11.3  FORMAT  MAP.) 

MAP 

46 

6 

FORMAT (lbHOOISPLAY  METHOD  I4,33H  IS  NOT  AVAILABLE.  USEO  METHOD 

1.) MAP 

47 

7 

FORMAT  C//15X, 2 8HTHE  OUTPUT  PRESENTATION  IS 

A) 

MAP 

48 

* 

FORMAT  <//15X, 25HTHE  CUANTITY  PRESENTED  IS) 

MAP 

49 

9 

F0°MAT(15X,43HA  COUNT  OF  CONTRIBUTING  DEPOSIT  INCREMENTS.) 

MAP 

50 

10 

FORMAT ( 15X v 42HEXP0SLRE  RATE  NORMALIZED  TO 

TIME  H*1  HOUR.) 

MAP 

51 

11 

FORMAT (15X , 24HEXP0SURE  PATE  AT  TIME  H+F10. 

1,9H  SECONDS.) 

MAP 

52 

12 

FORMAT <15X ,36HEXP0SL'RE  ACCUMULATED  BETWEEN 

TIME  H+F10 

.1, 22H  SECONDMAP 

53 

IS  ANO  INFINITY.) 

MAP 

54 

13 

FORMAT (15X , 36HEXPOSLRE  ACCUMULATED  BETWEEN 

TIME  H*F10 

•1,12H  AND 

TIMAP 

65 

1ME  H«-F10.1,9H  SECONCS.) 

MAP 

56 

14 

FORMAT <15X,60HTOTAL  MASS  ®ER  UNIT  AREA  OF 

CONTRIBUTING  OEPOSIT 

INCMAP 

57 

1REMENTS.) 

MAP 

58 

AH 


o  o  o  o 


15  FORMAT (15X,43HMASS  PER  UNIT  AREA  OEPOSITEO  BETWEEN  TIMES  F10.1,5H  MAP 


1ANO  F10.i»9H  SECONDS.)  MAP 

16  FOPMAT(/IX,F7.0,3X,2(10X,5H*****,F12.0,3X) ,20X,5H*****,/)  MAP 

17  FORMAT (15X , 41H ASSUMES  ALL  PARTICLES  ARE  GROUNOEO  BY  Tl.)  MAP 

15  F0SrtAT(15X,17HACTIV IT.Y  AT  TIME  F10.1,19H  OUE  TO  MASS  CHAIN  14)  MAP 

19  FORMAT (15X,26H MULTIPLE  BURST  BINARY  TAPE!  HAP 

20  FORMAT {15X.31HCR0UNr  ZERO  IS  LOCATED  AT  X  =  R10.1,5H  ,  Y  =  F10.1MAP 

1»  MAP 

?1  F0RMAT(1H1,41X,36HY-C00RDINATE  SCALES  FOR  SIDES  OF  MAP/1HQ)  MAP 

22  FO°MM  (//1X,F13»0»62X,F13.Q)  MAP 

23  FOR-  (15X , 46HTIME  (SECONDS)  OF  ONSET  OF  FALLOUT  DEPOSITION. »  MAP 


24  FORMA . <15X,50HTIME  (SECONDS)  OF  CESSATION  OF  FALLOUT  DEPOSITION.)  MAP 

25  FORMAT (15X» 50 HDIA ME  TER  (MICRONS)  OF  SMALLEST  DEPOSITED  PARTICLE.)  MAP 

26  F0OMAT«15X,49H0IAMETER  (MICRONS)  OF  LARGEST  DEPOSITED  PARTICLE.)  MAP 

27  FORMAT (l5Xf 58HMASS  CEPOSITEO  (<GM/H**2)  BY  PARTICLES  IN  THE  SIZE  RMAP 


1ANGE  ,E12.5,4H  TO  ,E12.5,  9H  MICRONS.)  MAP 

25  FORMAT (15X,77HH*1  HCUP  NORMALIZED  EXPOSURE  RATE  RESULTING  FROM  PARMAP 

1TICLES  IN  THE  SIZE  RANGE  ,E12.5,4H  TO  ,E12.5,9H  MICRONS.)  MAP 

29  FORMAT  ( 15X  , 25 HUN IT S  APE  ROENTGENS  PER  HOUR)  MAP 

30  FORMA T (15X , 19H UNITS  APE  ROENTGENS)  MAP 

31  FORMAT f 15X , 19 H UNITS  ARE  KGM/M**21  MAP 

32  FORMAT (15X »21HUNITS  ARE  CURIES/M**2)  MAP 

MAP 

?**»***»»**•**»»»»*»»»»******** MAP 

MAP 

99  IF(MAPRUN)  101,100, 101  MAP 

100  DO  1000  1=2,20  MAP 

FMTEXP(I)=BLAN<  MAP 

1000  FMTRUT(I)=BLA«i<  MAP 

TINC=5.0*OGX  MAP 

XCOORO=XMIN»OGX  MAP 

VINC= INC  MAP 

XCINC=VINC*OGX  MAP 

«L=1  MAP 

NX=NXMAP  MAP 

C  LEFT  IS  USED  HERE  AS  A  TEMPORARY  STORAGE  MAP 

LEFT= (XMAX— XI) /OGX  MAP 

C  PRINT  MAP  TITLE  MAP 

HRITE  (ISOUT, 7)  MAP 

C  SELECT  APPROPRIATE  CISPLAY  OPTION  COOE  MAP 

IF(JC(1))147,147,131  MAP 

131  IF (JC(1) -6 *132,132,147  MAP 

130  JC(1)=1  MAF 

132  N1=JC(1)  MAP 

GO  TO  (141,142,143,144,145,146) ,N1  MAP 

141  ASSIGN  150  TO  N2  MAP 

WRTTE  (ISOUT, 3)  HAP 

GO  TO  102  MAP 

142  ASSIGN  151  TO  N2  MAP 

WRITE  (ISOUT, 51  MAP 

GO  TO  102  MAP 

143  WRITE ( ISOUT ,19*  MAP 

ASSIGN  301  TO  N2  MAP 

IF(LREW.NE.O)  GO  TO  1431  MAP 

LREW=1  MAP 

REWIND  MBTAPE  MAP 

1431  WRITE  (MBTAPE)BITLUY  MAP 


59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 
76 
70 
BO 
51 
82 

83 

84 

85 
56 
87 
89 

89 

90 

91 

92 

93 

94 

95 

96 

97 

r»  • 

9C 

100 

101 

102 

103 

104 

105 

106 
107 
106 

109 

110 
111 
112 

113 

114 

115 

116 


58 


WRITE (M8TAPE) XMIN, XCAX , YMIN, YMAX, OGX, OGY 

MAP 

117 

Gn  TO  102 

MAP 

118 

c 

MAP 

119 

c**»* 

120 

144 

CONTINUE 

MAP 

121 

145 

CONTINUE 

MAP 

122 

146 

CONTINUE 

MAP 

123 

£♦*** 

124 

c 

MAP 

125 

147 

WRITE  (ISOUT, 6>N1 

MAP 

126 

GO  TO  130 

MAP 

127 

101 

KKL  =  1 

MAP 

128 

NX=NXMAP 

MAP 

129 

C 

LEFT  IS  USEO  HERE  AS  A  TEMPORARY  STORAGE 

MAF 

130 

LEFT* (XMAX-X11 /DGX 

MAP 

131 

GO  TO  1702 

MAP 

13? 

C  1C2 

PRINT  ORDINATE  DESCRIPTION 

MAP 

133 

C 

MAP 

134 

10? 

WRITE  (ISOUT,*) 

MAP 

135 

GO  TO  (161,162,163,164,165,166,167,169,169, 

17 1,172, 173, 174, 175, 176MAP 

136 

1,177,178,179,170,170) ,NPEO 

MAP 

137 

161 

WRITE  (ISOUT,9l 

MAP 

138 

GO  TO  170 

MAF 

139 

162 

WRITE  ( TSOUT, 101 

MAP 

140 

WRITE  (TSOUT, 291 

MAP 

141 

GO  TO  170 

MAP 

142 

163 

WRITE  (ISOUT, 111 T1 

MAP 

143 

WRITE  (ISOUT, 29) 

MAP 

144 

GO  TO  170 

MAP 

145 

164 

WRITE  (ISOUT, 12) T1 

MAP 

146 

WRITE  (ISOUT, 30) 

MAP 

147 

GO  TO  1?0 

MAP 

148 

165 

WRITE  (ISOUT, 13)Tl,T2 

MAP 

149 

WRITE  (ISOUT, 301 

MAP 

150 

GO  TO  170 

MAP 

151 

166 

WRITE  (ISOUT, 14) 

MAP 

152 

WRITE  (ISOUT, 31) 

MAP 

153 

GO  TO  170 

MAP 

154 

167 

WRITE  (ISOUT, 151T1.T2 

MAP 

155 

WRITE  (ISOUT, 31) 

MAP 

156 

GO  TO  170 

MAP 

157 

168 

WRITE  (ISO UT, 13)71, T2 

MAP 

158 

WRITC  (ISOUT, 30) 

MAP 

159 

WRITE (ISOUT ,17) 

MAP 

160 

GO  TO  170 

MAP 

161 

169 

WRITE  (ISOUT, 12*  T1 

MAP 

162 

WRITE  (TSOUT, 30) 

MAP 

163 

WRITE  (ISOUT, 17) 

MAP 

164 

GO  TO  170 

MAP 

165 

171 

WRITE  (ISOUT, l«)Tl,fASCHN 

MAP 

166 

WRITE  (ISOUT, 32) 

MAP 

167 

WRITE ( TSOUT ,17) 

MAP 

168 

GO  TO  170 

MAP 

169 

172 

WRITE  (ISOUT, 23) 

MAP 

170 

GO  TO  170 

MAP 

171 

173 

WRITE  (ISOUT, 24) 

MAP 

172 

GO  TO  170 

MAP 

173 

174 

WRITE  (ISOUT, 251 

MAP 

174 

59 

**Msw«u«wjau> 


o  o 


GO  TO  iro 

175  WRITE  (ISOUT,  261 
GO  TO  170 

176  write  cisour,27j  ti,t2 
GO  TO  170 

177  WRITE  (ISOUT, 281  T1,T2 
WRITE  (ISOUT, 29) 

GO  TO  170 


CODE  INSERTION  POINTS 


178 

170 


CONTINUE 

CONTINUE 


CODE  INSERTION  POINT! 


iro 

iao 


C 

c 


1701 

1702 
1021 
1022 


C 

C 


2023 


1023 

104 

103 

1031 


C 

C 


WRITE  ( ISOUT, 2  0 )  XGZ,YG7 
IF(JC(1).E0.3>  GO  TC  1702 

PRINT  a  PAIR  OF  PASTE-GN  V  SCALES  HERE 
WRITE  (ISOUT, 21) 

YY=YMlN*DGY*FLOAT(N>WAF! 

DO  1701  J=l, NYHAP 
WRITE  (ISOUT, 221  YY , YY 
YY= YY-OGY 

IF(LEFT-*VX1  1021,1022,1022 

NX=LEFT 

MM=NX/ (INC! 

jlj  j  ^  J 

LEFT  IS  USEO  HERE  AS  THE  NUM9ER  OF  PRINT  COLUMNS  IN  THE  LAST 
PRINTER  STRIP 
LEFT =NX- MM ♦(INC) 

IF  (LEFT.NE.91  GO  TC  2023 
M  =  MM 
LEFT  =  INC 
STRIPS 

00  110  ISTRIP=1,M 
MAPRUN=MAPPUN«-1 
IF  (JC (1) • EO. 31 
XC2=XCOORD*TINC 
XC3=XC2HINC 
WRITE  (ISOUT,  11HAPRLN 
WRITE  (ISOUT, 16)XC0CRD,XC2,XC3 
<1  * «L  M  N YM AP- 1 )  *  N  X  M A P 
IF(ISTRIP-MUO  3, 10  4,10  3 
KINC=LEFT-1 
VLEFT=LEFT 
XCIN=VLEFT»0GX 
GO  TO  1031 
KINC=INC-1 
XCIN=XCINC 
CONTINUE 
KLINK  =  KINCHl 

IF( JC (1) .EQ.3)  WRITE(MBTAPE) NYHAP, KLINK 
ROWS 

00  200  J=1,NYMAP 
KH=KL+KINC 
KOC=0 

DO  201  K=KL,KH 


GO  TO  1021 


MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
♦  MAP 
MAP 
MAP 
•MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 


175 

176 

177 
176 
17? 
160 
161 
18? 
193 

184 

185 
166 

167 

168 
1 39 

190 

191 

192 

193 

194 

195 

196 
107 
196 

199 

200 
201 
202 

203 

204 

205 

206 
20  7 
208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 
227 
226 

229 

230 

231 

232 


60 


f 

d- 


m 


m 


IF(OMAP(K»  .LT.CUTMAP»OMAP(KI=Q.O 

MAP 

233 

201 

FSUM=FSUM40MAP(K* 

MAP 

234 

C 

MAP 

235 

C 

NUM.  £RS  WITHIN  ROW* 

MAP 

236 

00  300  <=<L,<H 

MAP 

237 

Knc=Kocn 

MAP 

23  3 

C 

transfer  to  code  for  selected  presentation 

MAP 

239 

GO  TO  N2, (150, 151,301* 

MAP 

240 

c 

MAP 

241 

C  150 

CODE  FOR  POWER  OF  TEN  OISPLAY 

MAP 

242 

150 

IF  COMAP (IO  *105,106, 107 

MAP 

243 

105 

ASSIGN  121  TO  N3 

MAP 

244 

OMAP(K»=-OMAP(K* 

MAP 

245 

GO  TO  109 

MAP 

246 

107 

ASSIGN  300  TO  N3 

MSP 

247 

109 

H  =  ALOGIO ( OMAP( <) ) 

MAP 

248 

Hl=AMOD(H,l.Q* 

MAP 

249 

JMAP(K0C»=H-H1 

MAP 

250 

TF(JHAD(<r(C).rQ.O)  JMAP(<OC*=0 

Hup 

251 

FMTEXP  (KOCH*  =  FMT I 

MAP 

252 

FMTRUT (<DC*1I  =  FMT  F 

MAP 

253 

IF  ( JMAP (KDC* . NE. 0 ) GO  TO  1090 

MAP 

254 

JMAP(KCC»=0 

MAP 

255 

FMTEXP(KDC4l»  =  FMTA 

MAP 

256 

1090 

OMAP(K)  =  10.0**H1 

MAP 

257 

IF( OMAPf  <»- 9* 999)115, 115, 1091 

MAP 

258 

1091 

OMAP(K»=OMAP(<»  *10.0 

MAP 

259 

JMAP(<DC*=JMAP(<DC*  41 

MAP 

260 

FMTEXP(KOC4l*  =  FMTI 

MAP 

261 

GO  TO  115 

MAP 

262 

106 

JMAP(KOC»=0 

MAP 

263 

OMAP(K)=0.0 

MAP 

264 

FMTEXP(<DC41>  =  FMTA 

MAP 

265 

FMTPUT (KOC41*  =  FMTA 

MAP 

266 

GO  TO  300 

MAP 

267 

115 

GO  TO  N3, (300, 121* 

MAP 

268 

C  121 

PESET  SIGN  OF  MAP  COORDINATE 

MAP 

269 

1C1 

OMAPC  K)=-OMAP ( K» 

MAP 

270 

GO  TO  300 

MAP 

271 

C 

MAP 

272 

C  151 

CODE  FOp  TWO-LINE  FU.3  OISPLAY 

MAP 

273 

151 

JMAP  (<OC»=OMAP  (0/10.0 

MAP 

274 

ZMAP=JMAP(KDC» 

MAP 

275 

OMAP(<|=OMAP(<»-(ZMAP*10.0» 

MAP 

276 

FMTEXP (KOC  41*  =  FMTI 

MAP 

•*77 

FMTRUT (<OC4l)  =  FMTC 

MAP 

'■8 

FMTEXP(KOC4l»=FMTA 

MAP 

279 

FMTRUT (<OC4l*=FMTA 

MAP 

280 

300 

CONTINUE 

MAP 

231 

WRITEdSOUT ,2  )  (JHAP(<)  ,<=1,<DC» 

MAP 

282 

WRITE (IS0UT,4  >  (OMAPOO  ,K=XL,  KH» 

MAP 

283 

GO  TO  200 

MAP 

284 

301 

WRITE  (M8TAPE*  (OMAFOO  ,K=KL,  KH* 

HAP 

285 

200 

KL=KL-NXMAP 

MAP 

236 

IF  (JC(1*.EQ.3)  GO  TO  110 

MAP 

287 

WRITE  (ISOUT, 16JXCCCRO,XC2,XC3 

MAP 

2u3 

XCOORO=XCOORD4XCIN 

MAP 

289 

110 

KKL=KKL4lNC 

MAP 

290 

61 

«|.T^ 


— . 


oooooooooooooooooooooooooooooooooooooooon 


SUBROUTINE  PCHECK  PCHF<  1 

PCHEK  2 

THIS  SUP ROUTINE  DETERMINES  THE  TYPE  OF  MAP  REQUESTED  AND  PCHEK  7 

IT  INITIALIZES  FOR  THIS  HAP.  FOR  EACH  PARTICLc  IN  THE  DATA  BLOCK  PCHE<  A 

IT  CCHFUTES  THE  BCUAGRIES  OF  ITS  CGNTRI8UT10M  ELLIPSE  AND  PCHEK  5 

IT  LAPELS  IT  ACCORDING  TO  WHETHER  IT  HILL  CONTRIBUTE  TO  PCHEK  c 

SUBSEQUENT  MAP  COPE  LOADS  O'?  NOT.  IE  A  PAPTICLE  CONTRIBUTES  TO  PCHEK  7 

THE  CURRENT  HAP  COR’.’  LOAD  ,  SUBROUTINE  CALC  IS  CALLED.  PCHFK  8 

PC  HE !-'  B 

H.G.NORHENT  JUNF  23,1971  PCHEK  10 

PCHEK  11 

***i******#***»*********»****9»*»**  GLOSSARY  ********  ♦♦**#¥*****7*****»Pi~hEK  12 

PCHEK  13 

KTR (IP!  INDICATES  WHETHER  OR  NOT  THE  PARTICLE  IS  TO  »E  PCHEK  14 

CONSIDERED  IM  SUBSEQUENT  hftP  CORE  LOADS  -  -  PCHFK  IE 

0  -  CONSIDER  PARTICLE  SUBSEQUENTLY  PCHEK  16 

1  -  REJECT  PARTICLE  FOR  FURTHER  USE  PCHEK  17 

Y»P*«U  UPPER  Y  COORDINATE  LIMIT  FOR  PARTICLE  CONTRIBUTION  PCHEK  IP 


THIS  SUP ROUTINE  DETERMINES  THE  TYPE  OF  HAP  REQUESTED  AND 
IT  INITIALIZES  FOR  THIS  HAP.  FOR  EACH  PARTICLc  IN  THE  DATA  BLOCi 
IT  CCHFUTES  THE  BCUAGRIES  OF  ITS  CONTRIBUTION  ELLIPSE  AND 
IT  LABELS  IT  ACCORDING  TO  WHETHER  IT  WILL  CONTRIBUTE  TO 
SUBSEQUENT  HAP  COPE  LOADS  OR  NOT.  IE  A  PAPTICLE  CONTRIBUTES  TO 
THE  CURRENT  HAP  COR’.’  LOAD  ,  SUBROUTINE  CALC  IS  CALLED. 


H.G.NORMENT 


KTR (IP! 


JUNF  23,1971 


y°PMU 

XPRHU 

XPPNL 

YBRHL 

ASQ 


UPPER  X  COORDINATE  LIMIT  FOR  PARTICLE  CONTRIBUTION 

LOWER  X  COORDINATE  LIMIT  FOB  PARTICLE  CONTRIBUTION 

LOWER  Y  COORDINATE  LIMIT  FOR  PARTICLE  CONTRIBUTION 

SQUARE  CF  SEMI-AXIS  A  OF  THE  PARTICLE  CONTRIBUTION 

LIMIT  ELLIPSE 

SQUARE  OF  SEMI-AXIS  B  op  THE  PARTICLE  CONTRIBUTION 
LIMIT  ELLIPSE 

SIN  OF  THE  ORIENTATION  ANGLE  OF  THE  A  AXIS  OF 
THE  PARTICLE  CONTRIBUTION  LIMIT  ELLIPSE 
COSINE  OF  THE  ORIENTATION  ANGLE  OF  THE  A  AXIS  OF 
THE  PARTICLE  CONTRIBUTION  LIMIT  ELLIPSE 
LOG (BASE  El  OF  THE  RATIO  OF  THE  GAUSSIAN  PARTICLE 
CONTRIBUTION  DISTRIBUTION  MODE  VALUE  TO  QCUT 
COUNT  CF  AVAILABLE  PARTICLE  STORAGE  LOCATIONS  IN 
CORE.  THIS  IS  THE  NUMRER  OF  PARTICLES  REJECTED 
IN  PCHFCK. 

A  BLOCK  COUNT  OF  DATA  STORED  ON  TAPE  ANO/OR  IN  CORE 
MAGNITUOE  (I.E,  INTEGRATED  VALUE!  OF  A  PARTICLE 
PROPERTY  TO  BE  DISTRIBUTED  ON  THE  MAP 


PCHEK  14 
PC HF <  IE 
PCHEK  16 
PCHEK  17 
PCHEK  IP 
PCHEK  19 
PCHEK  20 
PCHEK  21 
PCHEK  22 
°CHpK  23 
PCHEK  ?4 
PCHEK  25 
PCHEK  26 
PPHEK  27 
PCHFK  25 
PCHEK  2® 
PCHEK  30 
PCHEK  31 
PCHEK  32 
PCHEK  33 
PCHEK  34 
PCHEK  35 
PCHEK  36 
PCHEK  37 
PCHEK  38 


ALSO  SEE  LINK®  GLOSS 

ARY 

PCHEK 

39 

PCHEK 

40 

►* PCHEK 

41 

PCHEK 

42 

COMMON  /SET17 

PCHEK 

43 

1CAY  yOETIDUE) 

,OIAMI201I 

tOMEAN 

yDNS 

,EXPO 

, PCHEK 

44 

20TTID (200) » IOISTR 

; IEXEC 

ylRISE 

yISIN 

yISOUT 

, PCHEK 

45 

3N0STR  ,  TID ( 20  0  I 

,SO 

ySSAM 

,THE 

yTMPl 

,  PCHEK 

46 

4TMP2  » T2M 

»u 

y  V  PR 

»W 

,HBURST 

y PCHEK 

47 

5SCL0H8  t NHOOO 

,ZV(20Q) 

yVX(200» 

,VYI200I 

PCHEK 

48 

COMMON  7PARDAT/ 

PCHEK 

49 

1X(500)  ,Y(500I 

,ZOUTT500> 

, SXOT (500) 

ySYOT (5001 

yRDUT (500) 

y  PCHEK 

50 

2PST500I  ,FMAS(500I 

,KTR(500) 

»  * 

y  GAMA 

,  BSQ 

,  PCHEK 

51 

3ASQ  »SIN A 

yCOSA 

y  WFMAS (2Q0) 

yYPRMU 

y YPRML 

, PCHEK 

52 

4TC500) 

PCHEK 

53 

COMMON  /RUNOAT/ 

PCHEK 

54 

INI J  ,NE 

»NREQ 

»NZ 

ylCTR 

,NXMAP 

y  PCHEK 

55 

2T1  ,T2 

t  HAPRUN 

,TGZ 

,IP 

yJC(16) 

y  PCHEK 

56 

3IC (18)  ,NYMAP 

,NTASK 

yNORO 

»XGZ 

,YGZ 

PCHEK 

57 

COMMON  /MAPOAT/ 

PCHEK 

58 

63 


o  O  OO  OO  OO  O  O  O  O  O  OO  OOt> 


,  XMIN 
,X2 


10MAP<15000) ,QCUT 
2XMAX 
3X1 

CCMMON/OECAY/ 

1  IGO 

2  , TEXIT 

common/output/ 

1  FTSNUM 


cCUTHAP 
>  YMAX 
,  MBTAPE 


» PGX 
,  Y  MIN 
,Z0EP 


,  jo 

,time 

.fp 


,KD0S 


C203>  ,FW 


,maschn 


,DGY 

,FSUH 


,  TENTER 


,  TT  A9 


♦DTLTAX 
, RUES AM 


,  JGO 


,PS!ZE  (2001  ,FMAS3{200>  ,PACTt200> 


210  FORMAT ( BHO  F0RMULAI6  ,67  H  IS  UNAVAILABLE.  COMPUTATION  HAS  CONTINUFO 
IFOR  A  REQUEST  OF  TYPE  6.1 

DATA  PR0GRM/6HPCHECK/ 

NE  =  0 

IF(IEXEC.EQ.l)  J  =  1 
IEXEC  =  0 

OO  7*7  IP=i,NIJ 

DETERMINE  IF  THE  OEFOSIT  INCREMENT  IS  GROUNDED 

IF((70UT(I°)-ZO£PI.L».1Q.OI  GOTO  75 
<T°< IPI=1 
GO  TO  777 

75  GO  TO <101, 102, 103, 104, 105,1 j6,107, 102,102.102,112,112,114 ,115,116 
1,117, 109, 110,  111, 12 Ct ,NREO 

101  COUNT  OF  GROUNDED  WAFERS 
101  F=FMASfIP» 

NORO=1 
GO  TO  100 

103  DOSE  RATE  AT  TIME  F*T1  SECONDS 

103  IF (T ( IP1  -111102,102,777 

104  OOSE  ACCUMULATED  FROM  TIME  H*Tl  SECONDS  TO  INFINITY 

104  IF(T(IP)-n»1041, 1041,1042 

1041  TENTEP=T1  -TGZ 
GO  TO  130 

1042  TENTER=TCIP>  -TGZ 
GO  TO  130 

105  DOSE  ACCUMULATED  FROM  TIME  H*T1  TO  TIME  H«-T2  SECONDS 

105  IF(T(IP)-T2»1051,777,777 

1051  IF«T(IP*-T1)1053, 1053, 1052 

1052  TENTER=T<IP!-TGZ 
GO  TO  130 

1053  TENTER=T1  -TGZ 
GO  TO  130 

106  TOTAL  PARTICLE  MASS  DEPOSITED 

106  F=FMASCIP> 

NORO=2 

GO  TO  100 


PCHEK  59 
PCHFK  60 
PCHF<  61 
PCHEK  62 
PCHEK  63 
PCHEK  64 
PCHEK  65 
PCHEK  66 
PCHFK  67 
PCHFK  66 
PCHF <  fig 
PCHEK  70 
PCHEK  71 
PCHEK  72 
PCHFK  7? 
PCHEK  74 
PCHFK  75 
PCHFK  76 
PCHEK  77 
PCHEK  7  3 
PCHEK  79 
PCHFK  3C 
PCHEK  81 
PCHFK  62 
PCHEK  63 
PCHEK  84 
PCHEK  35 
PCHEK  66 
PCHEK  37 
PCHEK  68 
PCHEK  89 
PCHEK  90 
PCHEK  91 
PCHEK  92 
PCHEK  93 
PCHEK  94 
PCHEK  95 
PCHEK  96 
PCHEK  97 
PCHEK  98 
PCHEK  99 
PCHEK100 
PCHEK101 
PCHEK102 
PCHEK103 
PCHEK104 
PC HEK10 5 
PCHEK106 
PCHEK107 
PCHEK106 
PCHEK109 
PCHEK110 
PCHEK111 
PCHEK112 
PCHEK113 
PCHEK114 
PCHEK115 
PCHEK116 


64 


s-tSl 


:  107  TOTAL  PARTICLE  HASS  DEPOSITED  BETWEEN  TIMES  T1  AND  T2  SECONDS 
107  IF(T<IP)-T2)l!I7l,777,777 

1071  TF(T<IP1-T1I777,777,106 

130  CALL  PAM2 

I  102  FIND  INDEX  CF  PARTICLE  SI7E  CLASS 
102  NOPO=2 

IF(N7.NE.0>  GO  TO  1020 

1022  IFtaestPStloi-PSIZEU))  .LT.  1.0E-2)  GO  TO  132 
J=  J+l 

IF ( J ,LE • IT AB)  GO  TO  1022 
CALL  FRPOR (PROGRM  ,-102,ISOUT» 

1020  JTA«J=TTABH 

DO  131  J=1 ♦ IT  AB 
X-JTAP-J 

IF  <**ACT  <tO  »GE.PS<IPH  GO  TO  1023 

131  CONTINUE 

CALL  ERROR  <PROGR*%  131,  ISOIJTI 
GO  TO  777 


IF (NREQ.EQ.4)  GO  TO  133 
F=FP  C  J  J*FMAS  (IP>  *’KFMAS(  J  1 
GO  TO  100 

F=FP(  j  >*FMAS1IP1 <FNASS(  J  J/RUFSAH 
GO  TO  100 


PCHF<117 
PCHEK118 
PCMEK11B 
°CHE<120 
PCHEX121 
PCHE<122 
DCHEX123 
PCHE<124 
OCHEK125 
PCHEK126 
POHFK127 
PCwn<123 
PCHE<12R 
PCHEK130 
PCHE<131 
PCHEK132 
PCHE<133 
PCHFK134 
PCHF<175 
PCHEK136 
PCHEK137 
PCHE<133 
PCHFK139 
PC.HF<140 
PCHFF141 
PCHE<142 
PCHEK143 
PCHE<144 
PCHEK145 
PCHEK14B 
PCHEK147 
PCHEK148 
PCHE<14B 
PCHE<150 
PCHEK151 
PCHEK152 
PCHEK153 
PCHEK154 
PCHEK155 
PCHE<156 
PCHEK157 
PCHEK158 
PCHEK159 
PCHEK160 
PCHEK161 
PCHEK162 
PCHEK163 
PC HEK164 
PCHEX1S5 
PCHEK166 
PCHE<167 
PCHEK168 
PCHEK169 
PC He <17  0 
PC  HE  KIT  1 
PCHE<172 
PCHEK173 
PC  HE KIT 4 


65 
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109 

110 

111 


CONTINUE 

CONTINUE 

CONTINUE 


C> 

C 


COOE  INSERTION  »0INTS 


120  CONTINUE 

WRITE  <ISCUT,210> NRfC 
NPF  0  =  6 
F=FMAS  < I Pi 


100  CONT' NUE 

COMPUTE  GAMA  AND  DETERMINE  THE  LIMITING  COORDINATES  OF  THE 
PARTICLE  CONTRIBUTION  ELLIPSE 

IF(F.LT.aCUT)  GO  TO  777 

GAMA  =  AL0G(  F/SXOT  (IP)/SYOT<IP)/QCUT/6,2e31<»531) 
IF(GAHA.LT.0.0*  GO  TO  200 
COSA=CCS(ROUT (IP) * 
sina=sin(rout  t ip* * 

AS0=  2.0*GAMA*SXOT t!P)**2 
BSO=  2.0*GAMA*SYOT tlP)**2 

YPPMU=Y(IPI  *  SOFT (ASG*SINA**2  ♦  RSG*C0SA«*2* 

YPPml  =  2.0*YCIP)-YPRMU 

DOES  THE  PARTICLE  CONTRIBUTE  TO  THE  MAP  WITHIN  ITS  VERTICAL 
«Y  AXIS*  LIMITS  - 


200 

205 


220 


230 

240 

777 


IF  (YPPMU.GT  ,Y1IN  ♦  CC-Y.  AND.  YPRHL.LT.  YMAX*  GO  TO  205 

<TR(IP»=1 

NE=ME*1 

GO  TO  777 

XPRHU=  X(IP)  #•  SORT  (ASQ*C0SA**2  *■  BSQ*SINA**2* 

DOES  THE  PARTICLE  CONTRIBUTION  LIE  COMPLETELY  BFYONO  THE  LEFT 
BOUNORY  OF  THIS  MAP  CORE  LOAD  - 

IF(XPRMU.LT.X1«-DGX)  GO  TO  200 
XPRML  =  2o  0  *XC IP)  -  XPRMU 

DOES  THE  PARTICLE  CONTRIBUTION  LIE  COMPLETELY  BEYOND  THE  RIGHT 
BOUNDRY  OF  THIS  MAP  CORE  LOAD  - 

IF { XPRML .LT »X2)  GO  TO  220 

KTPfIP*=0 

GO  TO  777 

WILL  THIS  G0NTRI6UTEP  ALSO  CONTRIBUTE  TO  SUBSEQUENT  MAP  CORE  L 

IFTXPRMU.GT.X2*  GO  TO  230 

KTRC IP*  =1 

NE=NE*1 

GO  TO  240 

KTRtIP*=0 

CALL  CALC 

CONTINUE 


PCHEK175 
Pr.HE<17f. 
PCHEK177 
•♦•’PCHEK178 
PCHEK17P 
PCHEK180 
PC HE<181 
PCHEK182 
PCHF<1«3 
PCHFK184 
•***PCHEK185 
PCHFK186 
PCHEK137 
PCHF<i88 
PCHEtfl89 
PCHE<190 
PCH£<191 
PCHFK192 
PCHFK193 

RCHF<194 

PrHE<195 

PCHPX196 

PCHEX197 

PCHEK199 

PCHEX199 

PCHEK200 

PC.HF<201 

PCHEK202 

PCHEX203 

PCHFK204 

PCHEK205 

PCHEK206 

PCHEK207 

PCHEK208 

PCHEK209 

PCHEK210 

PCHFK211 

PCHEK212 

PCHEK213 

PCHEK214 

PCHE<215 

PCHEK216 

PCHEK217 

PCHEK218 

PCHEK219 

PCHEK220 

PCHEK221 

PCHEK222 

PCHEK223 

0ADSPCHE<224 

PCHEK225 

PCHF<226 

PCHEK227 

PCHE<22« 

PCHEK229 

PCHFK230 

PCHEK231 

PCHEK232 
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SUBROUTINE 

POMP 

PDMP 

1 

POMP 

2 

THIS  SUBROUTINE  SORTS 

OUT  THOSE  PARTICLES 

THAT  HILL  CONTRIBUTE 

POMP 

T 

TO  SUBSEQUENT  HAP  CCRE  LOADS,  AND  DUMPS  THEM  ON  TO  T\ 

.'E  FOR 

PDMP 

4 

TEMPORARY 

STORAGE 

POMP 

5 

POMR 

f 

H.G.NORMENT  JUNE 

28,1971 

pnpp 

7 

POMP 

e 

♦  PDMP 

9 

priMP 

IP 

JL 

COUNT  CF 

FAPTIOLES  HOVEO  FROM 

UPPER  TO  LOWER  CORE 

POMP 

11 

(JL.LC.KP) 

POMP 

12 

JP 

COUNT  CF 

AVAILABLE  PARTICLE  STORAGE  LOCATIONS  »ASSEO 

i  POMP 

13 

IN  THE  PARTICLE  CORE  STORAGE 

BLOCK  SORT 

POMP 

14 

( JP.LE.NE.ANO. JP.LE.KP) 

POMP 

15 

<P 

NUMBER  OF  PARTICLES  IN  CORE  THAT  ARE  TO 

BE  DUMPED 

POMP 

16 

ONTO  TAPE 

PDMP 

17 

(KP=NIJ- 

NE) 

POMP 

18 

NE 

COUNT  CF 

AVAILABLE  PARTICLE  STORAGE  LOCATIONS  IN 

PDMP 

19 

CORE. 

THIS  IS  THF  NUMBER  OF 

PARTICLES 

REJECTEC 

POMP 

20 

IN  PCHECK. 

PDMP 

21 

NT  J 

A  BLOCK 

COUNT  OF  DATA  STORED 

ON  TAPE  AND/OR  IN  CORF 

POMP 

22 

POMP 

23 

ALSO  SEE  LINO  GLOSSARY 

POMP 

24 

ROMP 

25 

♦POMP 

26 

POMP 

27 

COMMON  /SET 1/ 

POMP 

28 

icay 

» OETIO  (12 ) 

, OIAM (2011  , OHEAN 

,9NS 

,EXPO 

,  PDMP 

29 

2DITID  (20D) 

.IOISTP 

, IEXEC  , TRISE 

,  ISIN 

,ISOUT 

,  POM® 

30 

3NOSTR 

i TIO (200) 

,SC  ,SSAM 

,TME 

,TMP1 

,  POMP 

31 

4TMP? 

*  T2M 

,U  ,VPR 

»W 

.HBURST 

,  POMP 

32 

SSCLOHB 

, NHOOO 

, ZV (2001  ,VX(200) 

,VY (200) 

POMP 

33 

COMMON  /PAOOAT/ 

POMP 

34 

1X(500) 

,Y(500) 

, ZOUT (500)  , SXOT (5001 

,SYOT (5001 

,ROUT (500) 

,  POMP 

35 

2PS (500) 

. FMAS (50  0 1 

,<TR(500)  ,  F 

,  GAMA 

,  RSQ 

,  POMP 

36 

3ASQ 

t  SIN  A 

,COSA  ,HFMAS(200) ,YRRMU 

,Y  RRML 

,  POMP 

37 

41(5001 

POMP 

38 

COMMON  /RUNOAT/ 

PDMP 

39 

INI  J 

»NE 

, NPEQ  , NZ 

,ICTR 

.NXHAP 

,  POMP 

40 

2T1 

t  T2 

,MAPRUN  ,TGZ 

,IP 

,JC  (13) 

,  PDMP 

41 

3IC (18) 

,  NY  MAP 

, NTASK  ,NORD 

,XGZ 

,YGZ 

POMP 

42 

COMMON  /MAPOAT/ 

POMP 

43 

10MAP(15000) ,CCUT 

,cutmap  ,ogx 

,DGY 

.DELTAX 

,  POMP 

44 

2XHAX 

,XNIN 

, YMAX  , YMlN 

,FSUM 

,RUFSAM 

.POMP 

45 

3X1 

»X2 

, MBTAPE  ,Z0EP 

PDMP 

46 

COMMON  /CONOAT/ 

POMP 

47 

1IPOUT 

» JPOUT 

.KPOUT  , <T  APE 

,ltape 

, MARRAY 

.PDMP 

48 

2NMAP 

.HXREO 

,IH  ,IV 

POMP 

49 

POMP 

50 

•POMP 

51 

POMP 

52 

OATA  PR0GRM/5HP0MP  / 

POMP 

53 

KP=NIJ-NE 

POMP 

54 

IF(NE.EO.O)  GO  TO  1000 

POMP 

55 

JP=0 

POMP 

56 

H=NIJ*1 

POMP 

57 

J=1 

POMP 

58 

68 


c 

c 

c 

c 

c 

c 


JL  =  0 


SORT  through  the  stcreo  particle  data  block  ano  hove  all 
PARTICLE  DATA  TO  «E  OUHPFO  INTO  LOWER  CORE  SO  THAT  IT  IS 
CONTAINED  IN  A  SOLIO  DATA  BLOC<  tl.E.  A  OATA  9L0C<  WITH  NO 
REJECTEO  PARTICLES  IN  IT» 

00  300  1=1, <P 
IF(KTRCI).EO.O)  GO  TO  300 
JP=JP*i 


POMP  59 
POMP  50 
POMP  61 
POMP  62 
POMP  63 
POMP  64 
POMP  65 
POMP  66 
POMP  67 
PnMP  66 


i 


00  200  K=J,NE 

POMP 

6° 

L=M-< 

POMP 

70 

IF{KTR1L).E0.1)G0  TC  100 

POMP 

71 

JL=JL+1 

POMP 

72 

K<=< 

POMP 

73 

C 

POMP 

74 

C 

MOVE  PAPTICLE  DATA  10  AVAILABLE  STORAGE  IN  LOWER  CORE 

POMP 

75 

C 

POMP 

76 

xm=xiL> 

POMP 

77 

y  <n=Ytu 

POMP 

76 

ZOUT(I)=ZOUT(D 

POMP 

7« 

r  m=T(D 

POMP 

90 

SX0T(I)=SX0T(LI 

POMP 

81 

SYOT  CI)  =  SYG'i‘(L1 

POMP 

92 

ROUT(T)=R'JUMLl 

POMP 

83 

PS(I»=PS«L1 

POMP 

94 

FMAS!l)=FMAS(L> 

POMP 

85 

GO  TO  260 

POMP 

96 

100 

JP=JP»1 

POMP 

67 

200 

CONTINUE 

POMP 

86 

250 

IPPOR=-250 

PDMP 

69 

GO  TO  2000 

PDMP 

90 

2b0 

J=K<*1 

PnMP 

91 

300 

CONTINUE 

POMP 

92 

IFIJP.LE.NEI 

GO  TO  400 

POMP 

93 

310 

IPOOR=-3io 

POMP 

94 

GO  TO  2000 

PDMP 

95 

400 

IFUP.LE.KP) 

GO  TO  FQO 

POMP 

96 

410 

IPROR=-Llu 

PDMP 

97 

GO  TO  2000 

POMP 

99 

500 

IF(JL.LE.<P)GO  TO  1000 

POMP 

99 

510 

IRR0R=-51Q 

POMP 

100 

2000 

CALL  ERROR  CPROGRM, IFROR,ISOUTI 

PDMP 

101 

1000 

WRITE! LT APES KP 

PDMP 

102 

WRITE  (LTAPE1CXII)  ,Y!I1  ,ZOUTm,T(I)  ,SXOT(I» 

,SYOT!I» ,ROUT  C II ,PS«I»,PDMP 

103 

lFMASdl  ,1=1, 

<P) 

POMP 

104 

RETURN 

POMP 

105 

END 

POMP 

106 

69 
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